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Abstract 

oo 

O ' The rare decay B — > Ji*(— !■ A' 7r)/i is regarded as one of the crucial channels for B 
physics as the polarization of the K* allows a precise angular reconstruction resulting in 
^ \ many observables that offer new important tests of the Standard Model and its extensions. 
^ \ These angular observables can be expressed in terms of CP-conserving and CP-violating 
quantities which we study in terms of the full form factors calculated from QCD sum rules 
on the light-cone, including QCD factorization corrections. We investigate all observables in 
the context of the Standard Model and various New Physics models, in particular the Littlest 
Higgs model with T-parity and various MSSM scenarios, identifying those observables with 
small to moderate dependence on hadronic quantities and large impact of New Physics. One 
important result of our studies is that new CP-violating phases will produce clean signals 
in CP-violating asymmetries. We also identify a number of correlations between various 
observables which will allow a clear distinction between different New Physics scenarios. 
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1 Introduction 



The penguin-induced flavour-changing neutral current (FCNC) transitions 6 — s and b ^ d 
are among the most valuable probes of flavour physics. They are characterized by their high 
sensitivity to New Physics (NP) contributions and the particularly large impact of short- 
distance QCD corrections to the relevant observables hke branching ratios and more local 
quantities, see Ref. [1] for a review. The decay 6 — 57 has probably been the most popular 
FCNC transition ever since its flrst experimental observation as 5 — )> K*'-/ at CLEO in 1993 
[2]. Despite its considerable success as benchmark probe, in connection with electroweak 
precision observables [3], its usefulness is limited by the number of observables it gives access 
to - the branching ratio and CP asymmetries like the time-dependent CP asymmetry in 
B K*-f [1]. 

Much more versatile in this respect is the decay b — s£~^i~ with the possibility to mea- 
sure, for instance, the differential decay rate in the leptons' invariant mass. One can also 
construct asymmetries, like the well-known forward-backward asymmetry {Aps), with differ- 
ing sensitivity to NP effects. A full angular analysis of B ^ K*{-^ K7T)£~^i~ would give 
access to a multitude of observables [5]. The downside of such measurements - low statistics 
- has started to be overcome at the B factories BaBar and Belle, with recent measurements 
of the forward-backward asymmetry in several bins in the lepton invariant mass and the K*^s 
polarization [BlEl- Current experimental results are compiled in Tab. [TJ The absolute num- 
ber of events observed is still rather small (230 at Belle [7]), making B — )■ K*i~^i~ one of 
the rarest B decays ever observed - at least if the resonance-dominated region around the 
charmonium resonances with B K*il){-^ is excluded. This situation will improve 

once the LHC experiments have started taking data, allowing one to probe the short-distance 
physics governing b — )■ si~^i~ at an unprecedented level of the angular spectrum: a recent 
study by the LHCb collaboration [10] predicts 7200 signal events (an improvement by an 
order of magnitude from the present situation) with a data set of 2fb~^, which corresponds 
to one nominal year of running. 

One difference between the experimental reach of B factories and LHC, though, is the 
preference of the latter for exclusive channels, mainly realized as B ^ K*£~^i~ . This implies 
that the analysis of this channel requires control not only over short- distance perturbative 
effects, described by Wilson coefficients in the relevant effective Hamiltonian, but also long- 
distance non-perturbative effects, described largely, but not completely, by form factors. It is 
the objective of our paper to provide such an analysis, based on QCD factorization [TT | [T2 | [T3] . 
including a full set of form factors calculated from QCD sum rules on the light-cone [H], 
and the dominant effects suppressed for large b quark mass. As the LHC has increased 
sensitivity to charged particles in the final state, we focus on the decays of neutral S's, 



Experiment 


BaBar [8] 


Belle [7] 


CDF [9J 


BR{B K*fx+fx-) X 10^ 


11.1 ± 1.9 ±0.7 


10.8+1:° ±0.9 


8.1 ±3.0 ± 1.0 


Number of BB events 


384 X 10^ 


657 X 10^ 





Table 1: Experimental results for the branching ratio of i? — t- i^*/i+/i~; the region around the 
charm resonances with B — )■ K*ij){-^ /i+/i~) is excluded. The first error is statistics, the second 
systematics. 
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^0 _^ K*0(^^ A'-7r+)/i+/i- and its CP-conjugate 5° K+n^)fi+fi- , which have the 

additional advantage that the flavour of the decaying B meson {B^ or B^) is unambiguously 
tagged by the flnal state. We also focus on i = fi which can be cleanly measured at the LHC; 
see Ref. [I5] for a discussion of /i O e effects. 

B — )■ K*n'^n~ decays have been investigated by many authors of whom we can cite only a 
few. In 1999, Ali et al. calculated the dilepton mass spectrum and Ap^ in the SM and various 
SUSY scenarios using naive factorization and QCD sum rules on the light cone [16]. Later 
it was shown by Beneke et al. [T2l [T3] that B — t- K*fi~^fi~ admits a systematic theoretical 
description using QCD factorization in the heavy quark limit rrih — )■ oo. This limit is relevant 
for small invariant lepton masses and reduces the number of independent form factors from 7 
to 2. Spectator effects, neglected in naive factorization, also become calculable. The drawback 
is that corrections to that limit are only partially known: for instance, in Ref. [T7] power- 
suppressed effects relevant for isospin asymmetries were calculated^ In Ref. [19], a calculation 
of -B — 7- K*fi'^fi~ using soft-coUinear theory (SCET) was presented. More recently, two 
analyses appeared which also use QCD factorization and focus on possible NP effects in CP 
asymmetries [20] and on observables available from angular distributions [21], respectively. 
There is also vast literature on NP analyses, with varying degrees of reliability of theoretical 
input for long-distance QCD effects and scope of observables considered. One intrinsic NP 
contribution, for instance, comes from an extended scalar sector. Most studies available so 
far, with the notable exception of Ref. [22], focus on the effects of such contributions on Aps, 
which turn out to be small, while we shall argue that the effect can best be seen in one 
particular angular observable not considered before, see Sec. [31 

In the present paper, we aim to improve on previous studies in the following way: 

• we include the full set of 7 form factors, rather than the 2 form factors in the heavy 
quark limit, calculated from QCD sum rules on the light-cone; we show that our set of 
form factors fulfllls all correlations required in the heavy quark limit, which has never 
been demonstrated before for any form factor calculation; 

• we give an up-to-date prediction of the B — )■ K*{^ K7T)fi~^fi~ observables in the SM 
and shall argue that the bulk of power-suppressed corrections is due to the difference 
between the full QCD form factors and their heavy quark limit; 

• we study all angular observables in the decay B K*{-^ Kn)^^ ^~ and identify those 
with small sensitivity to hadronic and large sensitivity to NP effects; 

• we include the effects of scalar and pseudoscalar operators, which are extremely sup- 
pressed in the SM, on all angular observables; 

• we study the effects of various NP models, including several manifestations of the MSSM 
and the Littlest Higgs model with T-parity. 

Our main results are collected in Sec. [71 

Our paper is organized as follows: in Sec. [2lwe review the theoretical framework, based on 
the trinity of effective Hamiltonian, form factors and QCD factorization. In Sec. [31 we discuss 
the (rather involved) kinematics of the decay and deflne the basic observables in the process. 
Sec. [H gives a short overview over the NP models whose effects we study. In Sec. [3 we deflne 

^ Very recently, BaBar has reported a positive result [18] for a deviation of the isospin asymmetry from the 
SM prediction, which so far, however, has not been confirmed by Belle [7]. 
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observables satisfying the requirements of theoretical cleanhness and high sensitivity to NP 
effects. Section [6l the centre part of our paper, contains the phenomenological analysis of 
those observables in the SM, in a model-independent way and in several selected NP scenarios. 
We conclude in Sec. [3 In the appendices we review the kinematics of four-body decays and 
show that for large b quark mass the form factors calculated, in Sec. [21 from QCD sum rules 
on the light-cone fulfill the relations imposed by heavy-quark symmetry. 

2 Theoretical Framework 

The theoretical framework which allows one to calculate the decay amplitude of B ^ K*fi~^n~ 
is quite involved and requires three different steps which are described in this section: 

• the separation of short- distance (QCD, weak interaction and new physics) effects from 
long-distance QCD in an effective Hamiltonian Hcs', 

• the calculation of matrix elements of local quark bilinear operators J of type {K*\J\B) 
(form factors); 

• the calculation of effects of 4-quark operators in Hcs which give rise to so-called non- 
factorizable corrections and can be calculated using QCD factorization (QCDF). 

QCDF is only valid for small invariant dilepton mass ~ 0(1 GeV^), or, equivalently, large 
K* energy E ~ 0{mB/2), which implies certain cuts on or E. In this paper, we restrict 
ourselves to 1 GeV^ < < 6 GeV\ The reasons will be discussed in Sec. 12.41 Obviously, 
all the above steps need to be under good control for a reliable prediction of the decay. We 
will discuss them in turn and also explain our strategy for calculating the B — )■ K*fi'^fi~ 
amplitude. 

2.1 Effective Hamiltonian 

The effective Hamiltonian for b — )■ sfi'^fi~ transitions is given by [23| El] 

^eff = (^t'^cff + -^""^cff ) (2-1) 

with the CKM combination Aj = VibV*. and 

6 

i=3 1=7,8,9, 10,P,S 

U^^ = C,{Ol - OX) + c^ipi - oi) . 

Although the contribution of "H^g is doubly Cabibbo-suppressed with respect to that of T^gg 
and hence often dropped, it proves relevant for certain observables sensitive to complex phases 
of decay amplitudes, so we keep it. The operators (9i<6 are identical to the Pi given in Ref. [23] . 
while the remaining ones are given by 

Or = -^m,{sa^,Pnb)F^\ O', = ^m,{sa,M)F'^'' , (2.2) 
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Os = -m,{m^,T'^PRh)G>'^\ O'^ = -m,{sa,,T-PLb)G^'''', (2.3) 

9 9 



O, = -(.-7^P^&)(/i7», = {sj^P^b)mrf^), (2.4) 



n \" fj,^ 1^" / \r~ r~/i o' 



e\ _ 



C^io = ^(s7^PL&)(/i7^5/i), Olo = ^(s7MPi?6)(/i7'^75/i), (2.5) 
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= Y^m,(sP^6)(/i/i), = —^m,{sPLb){fifi), (2.6) 

= T^'^b{sPRb){iJ.'y5iJ,), O'p = -—rmb{sPLb){iJ,'y5iJ,), (2.7) 

where g is the strong couphng constant and Pl,r = (1 T 75)/2. denotes the running 6 
quark mass in the MS scheme. The primed operators with opposite chirahty to the unprimed 
ones vanish or are highly suppressed in the SM, as are Os^p- We neglect the contributions 
of (9- for 1 < i < 6. These operators are generated in some NP scenarios, for instance in 
left-right symmetric models or through gluino contributions in a general MSSM, but their 
impact is either heavily constrained or turns out to be very small generically. 

The Wilson coefficients Ci in (12. ip encode short-distance physics and possible NP effects. 
They are calculated at the matching scale fi = mw, in a perturbative expansion in powers 
of as{'m^r), and are then evolved down to scales fi ~ mj, according to the solution of the 
renormalization group equations. Any NP contributions enter through Ci{mw), while the 
evolution to lower scales is determined by the SM. The inclusion of the factors IGn'^/g'^ = 
4:71 /as in the definition of the operators Oi>7 and the corresponding primed operators serves 
to allow a more transparent organization of the expansion of their Wilson coefficients in 
perturbation theory: all Ci are expanded as 



C. = Cr + ^cW+(^)Vf+0(aJ). (2.8) 



where C,-'^'* is the tree- level contribution, which vanishes for all operators but 02- In our 
normalization of operators also Cg°^ is non-zero. C]'"'' denotes an n-loop contribution. In our 
paper we aim at next-to-next-to-leading logarithmic (NNLL) accuracy, which requires the 
calculation of the matching conditions at yU = myy to two-loop accuracy. This has been done 
in Ref. [23]. NP contributions, on the other hand, will be included to one- loop accuracy onlyj§ 
Two-loop accuracy in the matching requires the inclusion of anomalous dimensions in the 
renormalization-group equations to three- loop accuracy. The corresponding O(a^) entries 
in the 10 x 10 SM anomalous dimension matrix have been calculated in Refs. [251 126]. On 
the other hand, the operators Og p are given in terms of conserved currents, i.e. they carry 
no scale-dependence, they do not mix with other operators and their Wilson coefficients are 
given by the coefficients at the matching scale. Og is also given by conserved currents, but 
mixes with Oi,...,6, via diagrams with a virtual photon decaying into Additional scale 

dependence in Cg comes from the factor 1/g'^. The latter dependence is also present in Cio, 
which otherwise would be scale independent. 



^An explicit calculation of two- loop corrections in the MSSM [24] shows that they are small. 
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In Tab. [2] we give all the SM values of the Wilson coefficients to NNLL accuracy. As we 
shall see below, in Eq. (13. ip . C^^g always appear in a particular combination with other Cj 

in matrix elements. It hence proves convenient to define effective coefficients Cj^g^ , and also 

C's icf ) which are given by [29] 



Att 1 4 20 80 



Cf^ — — Cs + C3 — - C4 + 2OC5 — ttCq , 
as b 6 

Cf = -Cg + Y{q'), 
as 

•-"10 ~ '-^10; '-"7,8,9,10" '-"7,8,9,10; {^■■~') 

as as ' ' ' 

with r = h{q^, rUc) (^^ Ci + C2 + 6C3 + 60Cj 

1/4 64 

- 2 ^(?'' ^b) [7Cs + 3 + 76C5 + y Q 

1/4 64 
--/^(g^0) C3 + -C4 + 16C5 + -Ce 



The function 



4 64 64 

3^^3 + yC^5 + ^C5. (2.10) 



arctan — . z > 1 



h{q ,mg) = -- In — 



9 \ u? 3 / 9^ ' 1,1 + VI - 2 

In ^ 2; < 1 

2 - 

(2.11) 

with z = 4mq/g^, is related to the basic fermion loop. 

We shall see below that B — A'*(— )■ Kir) fi'^ fi~ does not allow access to all the above 
coefficients separately: for instance, only the combinations Cs — C'g and Cp — C'p enter the 
decay amplitude. 



2.2 Form Factors 

The B — > K* matrix elements of the operators 10 5 p '^^^ t)e expressed in terms of seven 
form factors which depend on the momentum transfer q"^ between the B and the K* (g^ = 

{K*{k)\-s^^{l - ^,)b\B{p)) = -ze;{mB + mK*)Ai{q') + K^P " q),{^* " ?) ^'^^'^ 



ms + rriK* 
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n (,,\ 




(,,\ 


(,,\ 






Cj [fx) 


Cg w 






-0.257 


1.009 


-0.005 


-0.078 


0.000 


0.001 


-0.304 


-0.167 


4.211 


-4.103 


C'i(/i) 






C'4(/i) 




CM 








-0.128 


1.052 


0.011 


-0.032 


0.009 


-0.037 


-0.006 


-0.003 



Table 2: SM Wilson coefficients at the scale fi = nif, = 4.8 GeV, to NNLL accuracy. All other 
Wilson coefficients are heavily suppressed in the SM. The "barred" Ci are related to Cj as defined 
in Ref. [12]. Input: as{m]v) = 0.120, as(mb) = 0.214, obtained from as{mz) = 0.1176 [27], using 
three-loop evolution. We also use mt{mt) = 162.3 GeV [28], mw = 80.4 GeV and siv? 6w = 0.23. 

+ iq,{e* ■ q) ^ [A,{q^) - Ao(g^)] + e^,,.6*Vr ^^^^'^ , (2.12) 
q^ -' niB + rriK' 

with A,{q') = H^E^l^ A,{q') - A,{q') and Ao{0) = ^3(0); (2.13) 



{IC{k)\-sa^,q''{l+^,)h\B{p)) = ze^,,.e*Vr 2ri(g2) 
+ T,{q') [e;(m| - ml,) - (e* ■ q) {2p - q),] + T,{q'){e* ■ q) 



g2 



(2.14) 



with Ti(0) = T2(0). is the polarization vector of the K*. The form factors Ai and V are 
observables, i.e. scale independent, while the Tj depend on the renormalization scale /x. 
Aq is also the form factor of the pseudoscalar current: 

{K*\d^A^'\B) = {m, + m,){K*\si-i^h\B) = 2mK^{e* ■ q)A^{q^). (2.15) 

The form factors are hadronic quantities and call for a non-perturbative calculation. No 
lattice calculation of a full set of form factors is available yet. As a recent result we quote a 
(quenched) value for Ti(0) relevant for B K*-f: Ti(0) = 0.24 ± 0.03;o;oi [30]. Preliminary 
results from an alternative lattice calculation of Ti(0) have been reported in Ref. [31]. At 
present, a more promising method for calculating form factors at large energies of the final- 
state meson (i.e. at small g^) is offered by QCD sum rules on the light-cone (LCSRs) (s. 
Ref. [13] for reviews). This method combines standard QCD sum rule techniques with the 
information on light-cone hadron distribution amplitudes (DAs) familiar from the theory 
of exclusive processes [32]. It has been applied to B ^ K* form factors in, for instance, 
Refs. [331 133] • The key idea is to consider a correlation function of the 6—7-5 current and 
a current with the quantum numbers of the B meson, sandwiched between the vacuum and 
the K*. For large (negative) virtualities of these currents, the correlation function is, in 
coordinate-space, dominated by light-like distances and can be expanded around the light- 
cone. In contrast to the short-distance expansion employed in conventional QCD sum rules a 
la Shifman/Vainshtein/Zakharov [35], where non-perturbative effects are encoded in vacuum 
expectation values of local operators with vacuum quantum numbers, the condensates, LCSRs 
rely on the factorization of the underlying correlation function into genuinely non-perturbative 
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and universal hadron DAs (f). The DAs are convoluted with pro cess- dependent amplitudes Th, 
which similarly to Wilson coefficients can be calculated in perturbation theory, schematically 

correlation function ~ ^ ® (2.16) 

n 

The sum runs over contributions with increasing twist, labelled by n, and (8) means integration 
over the longitudinal momenta of the partons described by 0^"'-'. We shall see below that 
contributions of non- leading twist are suppressed by increasing powers of m/^./m;,. The 
same correlation function can, on the other hand, be written as a dispersion-relation, in the 
virtuality of the current coupling to the B meson. Equating dispersion-representation and 
light-cone expansion, and separating the B meson contribution from that of higher one- and 
multi-particle states, one obtains a relation (QCD sum rule) for the form factor. 
For B — )■ K* form factors the relevant correlation function is 

I J d%e--'Py{K*{p)\TJ,{0)jUym oc n(g2) (2.17) 

with Jb = di'j^b and = S7^(l —75)6, scr^,^g'^(l + 75)^ or si'j^b. The factor of proportionality 
contains four-vectors with open indices and/or mass factors like (m^ -|- mx*) etc., which 
are irrelevant for dynamics. LCSRs for all 7 form factors except for Aq are available at 
0{as) accuracy for twist-2 and-3 and tree- level accuracy for twist-4 contributions [3l]. For 
this paper, we have also calculated the LCSR for Aq to the same accuracy. The correlation 
function n(g^), calculated for unphysical p^, can be written as dispersion-relation over its 
physical cut. Singling out the contribution of the B meson, one has, for the pseudoscalar 
current = si'j^b, 

rrib — 

where J'b is the leptonic decay constant of the B meson, 

fsml = m,{B\bt^,d\0) . (2.19) 

In the framework of LCSRs one does not use (12.181) as it stands, but performs a Borel-trans- 
formation, 

5^ = ^exp(-t/M^), (2.20) 

with the Borel-parameter M^; this transformation enhances the ground-state B meson contri- 
bution to the dispersion-representation of 11. The next step is to invoke quark-hadron duality 
to approximate the contributions of hadrons other than the ground-state B meson by the 
imaginary part of the light-cone expansion of 11, so that 

5n^^ = ^ Ao{q') e-^'"^' + - rdtlmU^'^it) eM't/M') . (2.21) 
nib J-yl TT 

Subtracting the integral from both sides, Eq. f l2.2ip becomes the LCSR for Aq. sq is the 
so-called continuum threshold, which separates the ground-state from the continuum contri- 
bution. As with standard QCD sum rules, the use of quark-hadron duality above sq and 
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the choice of sq itself introduce a certain mo del- dependence (or systematic error) in the final 
result for the form factor. 

As an explicit example for a LCSR, we quote the tree-level result for Ti(0) as given in 
Ref. ESI: 



nib 



0±(^) 
2u 



•Jut) 



^(u) 1 



2u 

u — ai 



u du Jo Jo 2ai 



^ ^8^0 -4) 

Aia) + Via) 



2 /•! 



m 



b Juo 



i-^ <J uu(j)±{u) + 2Il{u) + uHs{u) 
2 du 



u pu 



dai I da2—{ S{a) - S{a) + t[^^ (a) - T:f^ («) + ^ (a) - Tf ^ [a] 
Jo «3 



1 d'^ . ^ ^ 



^(4), 



(2.22) 



no 



/i./2i(n) + /i^i?2(n) 



nib 



Rh(u) 



(2.23) 



where Uo is given by m^/so- /!:* and are the decay constants of, respectively, longitudi- 
nally and transversely polarized K* mesons. g^^'"'\ II and are DAs and integrals 
thereof, as defined in Ref. [31]. A, V, S", S and t/^-* are three-particle DAs. The precise 
definitions of all these DAs as well as explicit parameterizations can be found in Refs. [3Z]. In 
a slight abuse of language, we shall call Ri the twist-2 contribution to the sum rule, R2 twist-3 
and -R3 twist-4. u is the longitudinal momentum fraction of the quark in a two-particle Fock 
state of the final-state vector meson, whereas 01,2,3, with = 1, are the momentum frac- 
tions of the partons in a three-particle state. The light-cone expansion is accurate up to terms 
of order {niK* /ntbY- Up-to-date results for these DAs can be found in Ref. [22] • Although we 
only write down the tree-level expression for the form factor, radiative corrections are known 
for Ri [33] and the two-particle contributions to R2 [31] , and will be included in the numerical 
analysis. All scale-dependent quantities are calculated at the (infra-red) factorization scale 
ix^p = m\ — ml. The form factor itself carries an ultra-violet scale dependence. As a default, 
we choose /i = for that ultra-violet scale. 

It is clearly visible from the above formula that the respective weight of various contri- 
butions is controlled by the parameter niK*/fnb; the next term in the light-cone expansion 
contains twist-3, -4 and -5 DAs and is of order [mx* /nib)^- Numerically, the expansion works 
very well, with the 0{nij^*/nil) terms contributing less than 5% to the LCSRs. Note that the 
expansion is in niK*/nib only for = 0. For > 0, the expansion parameter is, see App. |Bl 
nibnix* /{ml — g^) ^ m^* /{2E), E being the energy of the K*. Obviously, the smaller E 
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B parameters 


fB m 


m 




200(25) MeV 


0.51(12) GeV 


2.2 GeV 


K* parameters 


fl 


/i*(2GeV) 


a|'ll(2GeV) 


a^'"(2GeV) 


220(5) MeV 


163(8) MeV 


0.03(3) 


0.08(6) 


quark masses 


mb{mb) |10] 


mc(me) [10] 


mt{mt) [28] 


4.20(4) GeV 


1.30(2) GeV 


162.3(1.1) GeV 



Table 3: Numerical values of hadronic input parameters, are parameters of the twist-2 K* 

DAs and are taken from Ref. [37], from where we also take all higher-twist parameters not included 
in the table. 

(and the larger g^), the more relevant the higher-twist terms. For E — rriK*-, the light-cone 
expansion breaks down. 

The LCSR method sketched above does not rely on mi, being a large (or hard) scale - 
LCSRs also work very well for D meson decays, see Ref. [38]. In order to calculate the 
B —J- K*fi'^fi~ decay amplitude, however, knowing the form factors is not enough: there are 
additional terms which can be calculated using QCD factorization, see the next subsection. 

As for numerics, we collect the most important input parameters in Tab. [3l We evaluate 
the sum rules at = 8 GeV^ and choose sq such that the minimum in is at 8 GeV^. The 
resulting sq lie all between 33 GeV^ and 37GeV^. Note that the LCSRs return results for /b 
times the form factor, rather than the form factor itself. Hence the LCSR must be divided 
by fs which, as one can see from Tab. |3l comes with a rather large error. It is well known 
that the resulting prediction for, for instance, Ti(0) is on the high side compared with the 
experimental result that can be extracted from the branching ratio of i? — )■ K*'j, assuming 
the absence of New Physics [12]. For this reason we fix the value of J'b to reproduce the 
experimental value Ti^^{0) = 0.268 [H]. This corresponds to setting Jb = 0.186 GeV - well 
within the allowed range quoted in Tab. |3l As the observables we calculate are ratios, the 
normalization and the precise value of cancel in the end - which is why we neglect the 
residual experimental and theoretical error of Tf^^(O) and Jb of about 7% [H]. The resulting 
values of the form factors at = are given in Tab. ID With Jb fixed, the errors of the form 
factors become rather small and are below 20%. 

As for the g^-dependence, it follows directly from the sum rules. In Fig. [1] we plot the 
central values of all form factors as functions of q^. 

2.3 QCD Factorization 

In addition to terms proportional to the form factors, the B — )■ K*fi^fi~ amplitude also 
contains certain "non-factorizable" effects that do not correspond to form factors. They are 
related to matrix elements of the purely hadronic operators Oi to Oq and the chromomagnetic- 
dipole operator with additional (virtual) photon emission. These effects can, in the com- 
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^i(O) 


^2(0) 


ViO) 


0.333 ±0.033 


0.233 ±0.038 


0.190 ±0.039 


0.311 ±0.037 


Ti(0) 


^3(0) 


^11 (0) 




0.268 ± 0.045 


0.162 ±0.023 


0.118 ±0.008 


0.266 ±0.032 



Table 4: LCSR results for = 0. T2(0) = Ti(0). The scale-dependent form factors Tj and _|_ 
are evaluated at /i = 4.8 GeV. The soft form factors y are introduced in Sec. 12.31 The error is 
calculated from varying sq by ±2 GeV^, by ±2 GeV^ and all hadronic input parameters according 
to their uncertainties given in Tab. [31 except for /b, see text. 




2468 10 2468 10 

(GeV^) (GeV^) 

Figure 1: Form factors from LCSRs for central values of input parameters. Left: Solid curve: 
long dashes: Ai, short dashes: A2, dot-dashed curve: V. Right: Solid curve: Ti, long dashes: T2, 
short dashes: T3. 



bined heavy quark and large energy limit, be calculated using QCD factorization (QCDF) 
methods [HI HSj [13]. Here large energy means large energy of the K*, E 0{mB/2). E is 
related to g^, the dilepton mass, by 

2mBE = m| ± m^. - . (2.24) 

For the phenomenological analysis in later sections, we require E > 2.1 GeV, which corre- 
sponds to < 6GeV^, well below the charm threshold. We would like to stress here that 
QCDF does not work for large above the charm resonances - here the only theoretical 
prediction we have are the contributions to the B — )■ K*n'^n~ matrix element given in terms 
of the form factors, which is probably a reasonable approximation at the 10 to 20% level. 

In the heavy quark and large energy limit, the number of independent form factors reduces 
from 7 to 2 which correspond to the polarization of the K* (transversal or longitudinal) and 
are usually denoted by and ^y. Neglecting for the moment 0{as) corrections, one can 
define the ^'s as [13] 

Uq') = ^^^V{q'), (2.25) 

Cii(g ) = — 7^ — M<i) Mq)- (2-26) 

Zll/ Til B 
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At this point we would like to recall that the large energy limit is also familiar from the 
description of hard perturbative QCD processes a la Brodsky-Lepage [12]. For instance, the 
electromagnetic form factor of the vr can be factorized into a convolution of the DAs of the 
initial and final state pion and a perturbative hard scattering amplitude, schematically 

i^.(Q') ~0.®Tj,®0,. (2.27) 

As in the previous subsection, the (S) stands for integration of the longitudinal momenta of 
all hard partons in the process. Note that is due to hard-gluon exchange between the vr 
constituents and hence is 0{as)- In this "hard mechanism" only valence-quark configurations 
contribute and all quarks have large longitudinal momentum. The subtlety with heavy meson 
decays, however, is that the power-counting of various contributions in l/m^ differs from that 
of the vr EM form factor and similar processes, and that a second mechanism contributes at 
the same (leading) order in l/m^. This is the so-called soft or Feynman mechanism, where 
not all partons in the meson participate in the hard subprocess; it involves highly asymmetric 
configurations where the spectator quark stays soft, but the light quark produced in the 
weak decay has large energy. This mechanism is described by the ^ form factors, which, as 
a consequence, are also called soft form factors, and the factorization formula for B — )■ K* 
decay form factors reads, schematically: 

F(g2) = D^{E) + <Pb®Th® ^k- + 0{l/mb) , (2.28) 

where D = 1 + 0{as) includes hard corrections to the weak vertex and E is the energy of 
the K*. As made explicit by the last term on the right-hand side, the above formula is not 
exact, but will receive corrections (both soft and hard) which are suppressed by powers of rrib. 
These corrections are unknown to date. An additional complication is that the separation 
between the "soft" contributions included in the form factor ^ and the "hard" contributions 
in the convolution is not clear-cut and requires the definition of a factorization scheme [TT] . 
In Ref. [13], the factorization scheme is defined by fl2.25p and f l2.26p . absorbing the hard and 
hard-spectator corrections into the definition of the .^'s. Eqs. fl2.25p and fl2.26p are then valid 
to all orders in a^. The redundancy of 5 form factors in the heavy quark limit induces a 
number of relations between them. We discuss these relations, and test their validity at finite 
b quark mass, in App. [Bl 

Coming back to the non-factorizable corrections to -B — > K*fi^fi~ mentioned before, it 
turns out that they can be included in a factorization formula very similar to that for form 
factors: apart from overall factors and the Lorentz structure, the relevant terms in the decay 
amplitude can be written as [13] 

r}'^ = iaC^^ + <Pb® T« ® 0,,^* + 0(l/mfe), (2.29) 

with a =_L, II and i = u,t. Note that the Ca^ in the above formula are not Wilson coefficients. 
They do, however, contain both factorizable corrections related to the rewriting of the full 
QCD form factors by the .^'s using f l2.28p and non-factorizable corrections related to the 
matrix elements of hadronic operators with virtual photon emission. Also note that 07,9,10 
only contain a 2-quark operator and hence do not induce non-factorizable corrections. 

2.4 Our Strategy 

Based on the above discussion, our strategy for calculating B — )■ /i*(— )■ KiT)fi~^fi~ decays is 
the following: 
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• we predict observables in the dilepton mass range IGeV^ < g2 < eCeV^; 

• we include the main source of power-suppressed corrections by using the full QCD 
form factors in the naively factorized amplitude, and the ^ form factors in the QCDF 
corrections; 

• we concentrate on the prediction of observables which are independent of the absolute 
values of form factors, and only depend on their ratios; 

• for the error analysis, we employ the correlated errors between form factors, which follow 
from the light-cone sum rules; 

• we include new-physics effects in the Wilson coefficients C7,9,io,5,p, and their primed 
counterparts, but not in the other Cj. 

A few comments are in order. Obviously QCDF breaks down close to the charm reso- 
nances. Technically, this shows up as a threshold at = 4m^. In order to stay sufficiently 
below the threshold, we set g^^^^ — 6GeV^. On the other hand, for small close to the 
kinematical minimum, the decay amplitude is dominated by the photon pole and by just one 
Wilson coefficient, C^^. Hence, by probing the region of very small g^, one does not get any 
new information as compared to the well-studied radiative decay b — >■ s'j. In addition, as the 
photon virtuality is small, there could be (unknown) resonance contributions from p or other 
mesons. In order to avoid this region, we set g^j^ = 1 GeV^. 

As for power-suppressed corrections, we view QCDF as an expansion in two small param- 
eters, l/rrih and and restrict ourselves to the first order in these parameters. While 0{as) 
corrections are completely covered by QCDF, those in l/nib are not included. One source of 
such corrections are obviously the differences between the 7 full QCD form factors and the 2 
soft form factors C,- Since the light-cone sum rules allow a calculation of all 7 form factors, 
all these corrections are included in our calculation. An obvious question, though, is whether 
these are all a^/rri}, corrections or not. As we shall see in Sec. 13.41 one QCDF correction to 
naive factorization is weak annihilation, where the quarks in the B meson annihilate. This 
contribution is of leading order in l/rrib, with unknown power-suppressed corrections. For 
the decay processes we consider, this contribution comes with small Wilson coefficients, so 
the impact of l/rrib corrections is negligible. Another potential source of l/rrib corrections are 
such with formal ag/mb counting, but an end-point divergence in the convolution integral. 
Such divergent integrals were found in power-suppressed corrections relevant for isospin vio- 
lation [TTj. They signal the breakdown of QCDF and call for soft contributions to mend the 
divergence. Such soft contributions involve soft gluons and hence as is not to be evaluated at 
a hard scale, but becomes non-perturbative, thus rendering this contribution 0{l/mb) in our 
counting. While it is not known how to calculate such contributions in the context of QCDF, 
similar contributions do occur in light-cone sum rules and are described by three-particle dis- 
tribution amplitudes of type {0\qGs\K*) , with G the gluonic field-strength tensor. Although 
these contributions could, a priori, be large due to soft-gluon effects, it turns out that such 
three-particle couplings are numerically small |37] and that, as mentioned in Sec. 12.21 their 
contribution to form factors is negligible. Based on this, we do not expect any sizeable effects 
of such terms at 0(1 /nib) and conclude that the main source of power-suppressed corrections 
are those from form factors. 
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3 Differential Decay Distribution and Spin Amplitudes 



In this section we discuss the kinematics of the 4-body decay B — )■ K*{^ Kti)^^ ^~ , define 
the angular observables in the spectrum and derive exphcit formulas in terms of form factors 
and Wilson coefficients. 



3.1 Differential Decay Distribution 

The actual decay being observed in experiment is not B — >■ K*fj,~^fj,~, but B — )■ K*{^ 
Ktt)^'^ ^~ . As discussed in Ref. [5], the additional information provided by the angle between 
K and vr is sensitive to the polarization of the K* and thus provides an additional probe of 
the effective Hamiltonian. 

The matrix element of the effective Hamiltonian (12. ip for the decay B — )■ K*{-^ K7i)^'^fi~ 
can be written, in naive factorization, as 



in 



{Knls^iCfPL + C','''PR)h\B) 



- '^{KTi\-sia^-q,{CfPR + C'fPL)h\B) 
+ (Kvr|s7^(CitPL + C'^^PR)b\B){-^^^^,^^) 



(3.1) 



+ {K'k\-s{CsPr + C'sPL)h\B){^i^^)+{K^:\-s{CpPR + C'pPL)h\B){ji^,fi)y 

To express the B — )■ Kn matrix elements in terms of the B — )■ K* form factors discussed 
in Sec. 12. 2^ one assumes that the K* decays resonantljEl. Then, one can use a narrow-width 
approximation by making the following replacement in the squared K* propagator: 

In this way, the form factors are independent of the K*K-k coupling Qk^k-k 0113], because it 
cancels between the vertex factor and the width 

Tk' = %^mK*f3', (3.3) 
48tt 



where ^ 

/3 = — — [m^, + 171% + - 2{m\,m\ + mj^ml + m^.m^)] ^^'^ . (3.4) 



Writing the matrix elements in Sec. 12.21 as 

{K*ik)\J,\Bip)) = e*''A,„ (3.5) 

where A^^j, contains the B — )■ K* form factors, the corresponding B — )■ Kir matrix element 
can then be expressed as 

{K{h)7r{k2)\J,\B{p)) = -DK^k^) WA,^, (3.6) 



^For a study of off-resonance effects, see Ref. 
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where [5] 



\DK'{e)\' = g^K.^^Sik'-ml.) = -^Sik'-ml.), (3.7) 

2 2 

^ i^A^ - ~ k^' = k1 + k^, K>' = k>^-k^. (3.8) 

With an on-shell K*, the decay is completely described by four independent kinematical 
variables: the dilepton invariant mass squared and the three angles 6k* , Oi and (p as defined 
in App. |X1 Squaring the matrix element, summing over spins of the final state particles and 
making use of the kinematical identities sketched in App. lAj one obtains the full angular 
decay distribution of 5° K*^{-^ i^"7r+ )//+//": 

d^T 9 

I{q\ei,eK*A). (3.9) 



dq^ d cos 6i d cos 6k* dcf) 327r 
where 

/(g^ Ou Ok* A) = A' sin^ Ok* + A'cos^ Ok* + {I2 sin^ Ok* + /a cos^ Ok*) cos20i 
+ I3 sin^ Ok* sin^ 6'^ cos 20 + J4 sin 20k* sin 2^; cos 
+ Is sin 2^x* sin Oi cos 

+ (Jg sin^ Ok* + Iq cos^ Ok* ) cos ^/ + sin 2^^* sin 6'; sin 
+ Jg sin 20k* sin 2^; sin + /g sin^ Ok* sin^ sin 20 . (3.10) 
The corresponding expression for the CP-conjugated mode 5° — Ji*°(— K~^iT~)fi~^ fi~ is 

7^7^ . = i^Ii<1^0,,0K*,<j^). (3.11) 

ag"^ ct cos t/; a cos c//^* a0 32-7r 

The function I{q'^,0i,0K*,4>) is obtained from (13.101) by the replacements [5] 

/-(«) , /-(a) , _ i 

-'1,2,3,4,7 ^-'1,2,3,4,7' -'5,6,8,9 ^ -'5,6,8,9' I^O.iZJ 

where /f"^ equals /j*''^'* with all weak phases conjugated. The minus sign in (I3.12p is a result of 
our convention that, while Ok* is the angle between the K*^ and the K~ fiight direction or 
between the K*^ and the K^, respectively, the angle Oi is measured between the K*^ (K*^) and 
the lepton /i~ in both modes. Thus, a CP transformation interchanging lepton and antilepton 
leads to the transformations Oi ^ Oi — ir and — >■ —0, as can be seen from Eqs. (lA.ip and 
(lA.2p . This convention agrees with Refs. [201113], but is different from the convention used 
in some experimental publications [10], where Oi is defined as the angle between K*^ and /i"*" 
in the decay, but between K*^ and /i~ in the decay. 

The angular coefficients J^- , which are functions of only, are usually expressed in terms 
of K* transversity amplitudes. Since we want to explicitly keep lepton-mass effects and 
include also contributions from scalar and pseudoscalar operators, this step deserves a closer 
look. 
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3.2 Transversity Amplitudes 

To introduce the transversity amplitudes, consider for the moment the decay B — K*V*, 
with the B meson decaying to an on-shell K* and a virtual photon or Z boson (which can 
later decay into a lepton-antilepton pair). The amplitude for this process can be written as 

Mim,n){B K*V*) = e],^,{m)M^,e*y'.{n) (3.13) 

where ey,{n) is the polarization vector of the virtual gauge boson, which can be transverse 
(ri = ±), longitudinal [n = 0) or timelike {n = t). In the B meson rest frame, the four basis 
vectors can be written as Il3l Il6] 



e^.(±) = (0,l,T^,0)/A (3.14) 

e^.(0) = (-g„0,0,-go)/v^, (3.15) 

e^.(t) = (go,0,0,g,)/V?, (3.16) 

where q'^ = (go, 0, 0, g^) is the four-momentum vector of the gauge boson. They satisfy the 
orthonormality and completeness relations 

ey^.(n)ey.^(n') = (3-17) 

Y,4'i^yv*in')9nn'=9''', (3.18) 

n,n' 



5 5 5 



where n, n' = t, ±, and gnn' = diag( 

The K* , on the other hand, is on shell and thus has only three polarization states, e^*(m) 
with m = ±,0, which read in the B rest frame 



e^,(±) = (0,l,±z,0)/A (3.19) 

e^.(O) = (A;„0,0,A;o)Mx*, (3-20) 

where = {ko,0,0, k^) is the four- momentum vector of the K* (note that k^ = —Qz)- They 
satisfy the relations 

e^.(m)ei^.^(m') = (3.21) 

J2 e*M^)e^K*{^') 6mm' = -g"" + (3.22) 

m,m' ^ 

The helicity amplitudes Hq, and can now be projected out from M^;^ by contracting 
with the explicit polarization vectors in fl3.13p . 

Hm = Mim,m)iB^K*V*), m = 0,+,-. (3.23) 

Alternatively, one can work with the transversity amplitudes defined as [15] 

= (i7+iTi7_i)/V2, Ao = Ho. (3.24) 
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In contrast to the decay of B to two (on-shell) vector mesons, to which this formahsm 
can also be apphed, there is an additional transversity amplitude in the case of i? — )■ K*V* 
because the gauge boson is virtual, namely 

A = M(odB^K*V*), (3.25) 

which corresponds to a K* polarization vector which is longitudinal in the K* rest frame and 
a V* polarization vector which is timelike in the V* rest frameEl 

If we now consider the subsequent decay of the gauge boson into a lepton-antilepton pair, 
the amplitude becomes 

M{B ^ K*V*i^ fi+fi-)){m) ^ e]^,{m) M,, Y^e*y^n)e1^,{n') g^r^' mipPL,Rfi). (3.26) 

n,n' 

This amplitude can now be expressed in terms of six transversity amplitudes y q and y g, 
where L and R refer to the chirality of the leptonic current, as well as the seventh transversity 
amplitude At. The reason that for At no separate left-handed and right-handed parts have to 
be considered can be seen as follows. Noticing that the timelike polarization vector in fl3.16p 
is simply given by eyt,{t) = q^/ \/^, one can see from current conservation, 

g^(/i7» = 0, g^(/i7''75/i) = 2zm^(/i75/i), (3.27) 

that the timelike component of the V* can only couple to an axial- vector current. In addition, 
this shows that At vanishes in the limit of massless leptons. 

Now, having shown that the amplitude of the sequential decay B K*V*{^ l^^ 
can be expressed in terms of seven transversity amplitudes, it is clear that this is true for 
all contributions of the operators Og^ and 0{q to the decay of interest, B — )• i^*(— )■ 
K7i)fi'^fi~ , regardless of whether they originate from virtual gauge boson exchange (i.e. photon 
or Z penguin diagrams) or from box diagrams. 

Does this also apply to decays mediated not by a vector, but a scalar and pseudoscalar 
operator? Inspecting Eqs. (12.61) . (12.71) and (13.271) . one can see that the combination (Op — O'p) 
can be absorbed into the transversity amplitude At, because it couples to axial- vector currents, 
just like the timelike component of a virtual gauge boson. However, this is not possible for 
the scalar operators 0^^\ Therefore, the inclusion of scalar operators in the decay B — )■ K*{-^ 
K7r)fi'^IJ,~ requires the introduction of a an additional, "scalar" transversity amplitude, which 
we denote As- 

To summarize, the treatment of the decay B — > K*{-^ Kti)^^^~ by decomposition of 

L Ft 

the amplitude into seven transversity amplitudes A_^ y g and At is sufficient as long as the 

operators o!j \ and Cp are considered, but has to be supplemented by an additional, eighth 
transversity amplitude As once contributions from scalar operators are taken into account. 

Finally, we give the explicit form of the eight transversity amplitudes (up to corrections 
of 0{as), whose discussion we postpone until Sec. 13. 4p : 



A^L.R = nV2X'/^ 



thb -r rrix* q 



(3.28) 



^Unlike sometimes stated in the literature, At does not correspond to a timelike polarization of the K* 
meson. As mentioned above, the K* decays on the mass shell and thus has only three polarization states. 
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AiiL,B = -NV2{ml - ml,) 



[iCf - Cf) T (Cit - Ct)] 



iriB — rriK* 



^{Cf - Cf)T,{q') 



(3.29) 



N 



OL,R 



2mK* V 



[{Cf - Cf) T (Cit - Cf)] 



X 



"^1 - "^Ic" - q'^){rnB + mK*)Ai{q'^) - A 
m| + 3m^*-g2)T2(g2) 



2mb{Cf - Cf 



A2{q') 
rriB + rriK* 

A 



2 2 



Tsiq') 



A, 



N 



yl/2 



2{Cll-Cf) + ^{Cp-C',) 



A5 = -iVAV2(C5-C;,Mo(g2), 



where 



1/2 



(3.30) 

(3.31) 
(3.32) 

(3.33) 



with A = + m^, + — 2[m\m\* + m\-,q^ + m\q^) and = yl — Arn?^/q^. 
3.3 Angular Coefficients 

With the eight transversity amphtudes defined in the preceding subsection, the angular coef- 
ficients li in (13.1 op can be written as 



I ^f^u) r|^L|2^|^L|2^(^^^- 



^ O^il' + \A\? + {L^ R)] + ^Re [AiAf + A\Af) 



(3.34) 





Am 

\A'.? + + ^[\At\' + 2Re(A^O] + 


(3.35) 


^2 = 


f 0^ir + i4l' + (^^^)]' 


(3.36) 


^2 = 


-f3l[\A^o\' + iL^R)], 


(3.37) 


/3 = 


l^l[\Ai\'-\Al^\' + {L^R)], 


(3.38) 


h = 


-^(31 [Re(AoMf *) + (L ^ i?)] , 


(3.39) 
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(L R) 



II = 2/3^ [Re{A\Ai*) - (L -> R)] 



4/3, 



Re [A^A*s + (L ^ R)] 



m, 



Re{AiA*s + A!l'A*s 



R A*' 



(3.40) 

(3.41) 
(3.42) 



lm{A^AF) -{L^R) + ^ Im(A^A^ + A«A 



m, 



/s = --j=Pl [lm{A^Ai*) + {L-^ R)] , 
h = /3l [MAfAi) + (L R)] . 



A few comments are in order: 



(3.43) 

(3.44) 
(3.45) 



In contrast to the transversity amplitudes themselves, the angular coefficients /j are 
all physical observables. In fact, they contain the complete information that can be 
extracted from the measurement of the decay — )■ K*^{^ K"'K~^)fi~^fi" . We will 
discuss in Sec. [5] which combinations of the angular coefficients constitute theoretically 
clean observables. 



In the limit of massless leptons, the well-known relations = 3/| and 



-II hold. 



The coefficient Jg vanishes unless contributions from scalar operators and lepton mass 
effects are taken into account. Therefore, to our knowledge, it has never been considered 
in the literature before. However, it is a potentially good observable for scalar currents. 
We will come back to this point in Sec. 16.2.31 



3.4 Additional Corrections to Transversity Amplitudes 



As mentioned in Sec. 13.21 the transversity amplitudes f l3.28p to f l3.32p do not include effects 
from spectator interactions, which do induce, on the one hand, 0{as) corrections and, on the 
other hand, corrections from weak annihilation (WA). These corrections have been calculated 
within the QCD factorization (QCDF) framework in Refs. [12] and [13] in terms of the soft 
form factors C,± and discussed in Sec. 12.31 

In Ref. [12], there are two types of 0{as) corrections, factorizable and non-f actor izable. 
The factorizable corrections arise when expressing the full form factors in terms of and 
are given by the radiative corrections in Eqs. (IB.ip -f lB74|) and therefore are redundant in our 
set-up. The only exception arises upon expressing the running b quark mass in the operators 
Eqs. (12. 2p and (12. 3p . by a mass parameter in a different renormalization scheme. In 
the numerical analysis, however, we use the running b quark mass in the MS scheme, so all 
factorizable 0{as) corrections calculated in Refs. [121 US] have to be dropped. 

The second QCDF correction to the transversity amplitude in Sec. I3.2l is given by the WA 
contribution, T||°2(m,w) in the notation of Ref. [12]. It is induced by the penguin operators 
O3 and O4 and hence is numerically small, see Tab. [2J This is a term which is leading in l/mi, 
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and 0(as), so in principle one should also include power-suppressed and radiative corrections. 
However, in view of its small size, we feel justified in neglecting them. As discussed in Ref. |13j . 
there are further WA corrections which are suppressed by one power of rrib with respect to 
the leading terms. For the leading CKM amplitude in At, see Eq. (12. ip . these are again due 
to penguin-annihilation diagrams and hence can be neglected due to the smallness of the 
Wilson coefficients. The WA contribution to the amplitude vanishes for — )■ K*^ix^ ix~ , 
but contains, for 5+ — > K*^fi'^fi~, the factor C2 ~ 1 and hence should be included for this 
process. As, in this work, however, we focus on neutral B meson decays, we can neglect all 
WA contributions except for Ty'^^l^? i^)- 

On introducing the chirality-flipped operators, the T^*||"^ introduced in Sec. 12.31 are pro- 
moted to 7^11*'"'' corresponding to the notations of Ref. [H] . In terms of these quantities, we 
can define the additional corrections to the transversity amplitude^: 

AA?^ J = V2n'^-^ (ml - g2)(^+W,WA+nf ^ ^^^+(.)^ 
A^fg^ = -V2N ^ (ml - ,2)(r-(*)'^^+'^^ + X^T^^^^) , 

,2 ^2\2 



m^^.ml^ 



AA^gg^ = ^'^""^^ % m,(7^^^^^'^^™ + XuTf'''^) . (3.46) 



The superscript, WA+nf, on 7]^^*"* indicates that only contributions from WA and non- 
factorizable 0(«s) corrections are to be included. In accordance with Ref. [20], we define 
Au = Xu/Xt- The total transversity amplitudes are given by the expressions in (I3.28I) - (I3.30I) 
plus the above terms AA^*^^^ . Note there are no corrections to At or As- 



4 Testing the SM and its Extensions 
4.1 Preliminaries 

The multitude of observables accesssible in 5 — )■ K*fi'^fi~ decays allows one to test the 
SM and its extensions more locally than is possible through global quantities like branching 
ratios and the dimuon mass spectrum. The goal of this section is to describe very briefiy the 
extensions of the SM that we will analyse numerically in Sec. El To this end we distinguish 
different classes of models using two properties: 

• the presence or absence of additional operators in the effective weak Hamiltonian relative 
to the SM ones, 

• the presence or absence of new sources of flavour and CP violation beyond the CKM 
matrix. 

When appropriate, we will also comment on the correlation between observables in 5 — )■ 
K*fi~^fi~ and other important observables in B physics, such as the mass difference in the 

^It should be noted that the functions F^'^^f^ entering the non-factorizable corrections are defined with a 
different overall sign in Refs. [13] and [48] . 
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neutral B^^s meson systems, AM^^^, and the time-dependent CP asymmetries 5* in various 
decay channels. These are: S^Ks^ measured in — )■ J /ipKs, which in the SM equals sin 2/3, /3 
being one of the angles of the unitarity triangle; S^Ks^ originating from the b — )■ sss penguin- 
decay Bd — )■ (pKs, which in the SM also equals sin 2/3, but is sensitive to new CP-violating 
phases in 6 — )■ s transitions; 5*^^, measured in Bg — )■ J/'ip(p, given by the Bg mixing phase, 
which in the SM is close to zerojj 

First, however, we would like to stress the importance of the observables in i? — )■ )■ 
K7i)^~^fi~ for tests of the SM. 

4.2 Standard Model 

The importance of the observables discussed in the present paper for tests of the SM originates 
from the following facts: 

• Several of the observables we consider are predicted to be strongly suppressed in the 
SM or even vanish so that New Physics (NP) effects can be seen more easily than in 
the branching ratio of 5 — K*fi~^fi~ which is measured to be consistent with the SM 
expectations. 

• The relatively small number of relevant SM parameters which are already well con- 
strained by a number of processes allows rather definitive predictions for many observ- 
ables subject mainly to the theoretical uncertainty of form factors. 

• In certain cases the sign of a given observable has a unique prediction in the SM, which 
can be tested more easily than the magnitude itself. 

However, to use these facts in a meaningful way, it is essential to have reliable calculations of 
the relevant form factors. In fact the use of the improved form factors presented in our paper 
and the consideration of correlations between the uncertainties of the different form factors 
allows one to obtain rather reliable predictions for angular coefficients in the SM. 

4.3 Models with Minimal Flavour Violation (MFV) 

The simplest class of extensions of the SM are models with constrained MFV (CMFV) in 
which the operators are SM-like, all flavour violating transitions are governed by the CKM 
matrix and also CP- violating observables are SM-like [501 Ell [52] . NP effects in this class of 
models affect the Wilson coefficients at the scale 0{mw) and collect NP contributions from 
scales higher than mw- At the same time, the QCD renormalization group evolution down to 
scales lower than mw is universal for all CMFV models and the same as in the SM. It can be 
shown that this class of models is characterized by strong correlations between observables in 
Bd, Bs and K processes [51] . 

The implications of these correlations and of the fact that all flavour violating interactions 
are governed by the CKM matrix are striking: deviations from the SM predictions for most 
weak-decay observables are bounded to be at most 50%, often much less [S2l [M]- Moreover, 
the predictions for CP asymmetries are generically identical to the SM ones. Consequently, 
the distinction between various models of this class on the basis of global quantities like 

^There are, however, hints from the Tevatron that this phase might actually be large [15] . 
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branching ratios is very challenging^ We will investigate whether more local quantities like 
the angular observables considered in our paper could be helpful in this respect. Obviously 
any measurement of a non-standard CP-asymmetry would immediately signal new sources of 
CP violation beyond the CMFV framework. 

The symmetry-based definition of MFV [58] does not exclude the appearance of addi- 
tional non-SM operators. Compared to the CMFV framework, large effects can be expected 
in particular from scalar operators. The most popular model in this respect is the MFV 
MSSM which boasts an extensive literature. In this model the CKM matrix remains the 
only source of flavour and CP violation, but, in particular at large tan /3, additional scalar 
and/or pseudoscalar operators, not present in CMFV models, enter the game. For instance in 
Bs,d /^^/^~ decays the presence of such scalar or pseudoscalar operators which are induced 
by Higgs penguin diagrams can enhance the branching ratio by one order of magnitude with 
respect to the SM, models with CMFV and the Littlest Higgs model with T-parity. 

The discovery of Bg — )■ ^~ with a branching ratio of 0(10"*^) would be a clear signal 
of non-standard scalar or pseudoscalar operators. In Sec. E] we shall demonstrate that it is 
precisely the angular observables in B ^ K*{-^ K7r)fi'^fi~ that allow one to distinguish 
whether scalar or pseudoscalar currents would be responsible for such an enhancement. 

4.4 Flavour-Blind MSSM 

As a modest, but interesting modification of the MSSM with MFV we consider the so-called 
Flavour Blind MSSM (FBMSSM) |59l |60l iB [62]. In this framework, the CKM matrix 
remains the only source of flavour violation, while new flavour conserving, but CP-violating 
phases are present in the soft sector. To be more specific, we assume universal soft masses for 
the different squark flavours and flavour-diagonal trilinear couplings at the electroweak scale, 
allowing for complex values of the latter parameters. 

One would naively expect that CP-violating, but flavour conserving observables, such as 
electric dipole moments (EDMs), would be the best probes by far for CP violation in the 
FBMSSM. However, it turns out that large new CP-violating effects in flavour physics can 
still occur at an experimentally visible level. In particular, there is still a lot of room for CP- 
violating asymmetries in B decays, like Acp{b — )■ 57) and S^Ks^ while S^^ is well constrained. 

The dipole operators play the most relevant role in this context and, as shown in Ref . |62] , 
striking correlations between electric dipole moments of neutron and electron, the Sf^Ks asym- 
metry and Acp{b — )■ 57) are present in this model. In particular this framework links the 
explanation of the suppression of Sif^Kg relative to S^Ks^ measured at the B factories, to a 
large enhancement of Acp{b — )■ 57) and the EDMs of the neutron and the electron, (i„ g, over 
their SM values. 

Therefore it is of interest to investigate whether the addition of flavour conserving, but 
CP- violating phases in the soft sector would have a visible impact on B K*{-^ K7r)fi^fi~ 
and consequently whether the FBMSSM could also be tested with the help of the observables 
discussed in the present paper. As we will see below, very interesting and predictive results 
are in fact obtained in this framework. 

''However, in the case of S ^^7, B Xst^i~ [55l [56| and in particular in B ^ K*£+£~ [57] sizable 
effects have been identified in the ACD model with one flat universal extra dimension for low compactification 
scales. 
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4.5 Littlest Higgs Model with T-Parity (LHT) 

Another class of models of interest are those in which the operators remain as in the SM, 
but new sources of both flavour and CP violation beyond the CKM matrix are present. In 
this class of models the CMFV correlations between B^, Bg and K observables are generally 
violated and much larger NP effects than in CMFV models are possible. 

A prominent example of this class of models is the Littlest Higgs model with T-parity 
[63| IM] in which the interactions between SM quarks and heavy mirror quarks, mediated 
by new heavy charged and neutral gauge bosons, involve a new mixing matrix that differs 
from the CKM matrix [65] and is parameterized by three new mixing angles and three new 
CP- violating phases |66] . 

A number of detailed analyses of FCNC processes in the LHT model has shown that large 
departures from SM predictions for FCNC processes are still possible in this model while 
satisfying all existing constraints [671 EHl EHl [70]. In particular the CP asymmetry S^^ can 
be enhanced by an order of magnitude relative to the SM prediction [671 E] which would be 
welcome if the data from the Tevatron [19] will be confirmed by more accurate measurements 
at LHCb. We will investigate whether the LHT model can also be tested efficiently by means 
of the angular observables in i? — )■ K*{-^ K7r)fi~^fi~. All loop functions with mirror quarks 
and new heavy weak-boson exchanges have been calculated in Refs. [671 [68] . A very recent 
paper, Ref. [70], finds additional contributions to Z penguin relative to Ref. [M]- We will 
investigate the importance of these terms in our analysis. 

4.6 General MSSM 

Finally we also consider the MSSM with generic flavour- and CP- violating soft SUSY-breaking 
terms. In such a framework one is confronted with a large number of free parameters which 
make it very difficult to perform global analyses. 

The flavour-mixing off-diagonal entries in the squark mass matrices, usually called mass 
insertions, present in this framework can lead to complex contributions to the Wilson coeffi- 
cients of all operators in Eqs. (12. 2p to (12. 7p . On the other hand, the mass insertions are not 
completely free parameters, but are constrained by measurements of many FCNC processes 
like BR{B X,7), BR{B AM,, AM^, S^Ks and others [721 [731 [72 [Z5] • The 

remaining parameter space still allows sizeable effects in the Wilson coefficients governing 
B — )■ K*fi^fi~. The general MSSM contributions to these Wilson coefficients have been stud- 
ied for both inclusive decays, B — )■ Xsfj.~^fi~ [761 and exclusive channels, B — )■ K^*^fj,'^fj,~ 
[171 [78] . In Sec. 16.3.41 we investigate the possible impact of these contributions on the ob- 
servables discussed in the following Sec. [51 focusing in particular on the question of how to 
distinguish the general MSSM framework from the other models described above. 

5 Observables 

As discussed in Sec. [3l the decay B^ — )■ A'*°(— K^n^)fi~^fi~ is completely described in terms 
of twelve angular coefficient functions J^-'^''. The corresponding CP-conjugate mode B^ — )■ 
A'*°(— )■ K~^7i~)fi~^fi~ gives access to twelve additional observables, the CP-conjugate angular 
coefficient functions 1^""^ These quantities have a clear relation to both experiment and theory: 
theoretically they are expressed in terms of transversity amplitudes, and experimentally they 
describe the angular distribution. A physical interpretation of these ij^""^ can be drawn from 
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SM 


18 


22 


SM + 


20 


24 



Table 5: Number of independent observables in i? — t- )• Ktt)^^^ , depending on whether 

lepton mass effects and/or scalar operators are taken into account. 

Eqs. (I3.34P to f l3.45p . For example, Jg depends on scalar operators and J7 to Jg depend on the 
imaginary part of the transversity amplitudes, and consequently on their phases, which come 
either from QCD effects and enter the QCD factoriation expressions at ©(a^), see Sec. [21 or 
are CP-violating SM or NP phases. 

To separate CP-conserving and CP-violating NP effects, we find it more convenient to 
consider the twelve CP averaged angular coefficients 

= (/<°' + /<°') /'-^ (5,1) 

as well as the twelve CP asymmetrie^ 

^^{lf_jf)l^±Il, (5.2) 

These are our primary observables that will be used in the phenomenological analysis in Sec. El 
They offer a clean and comprehensive way to analyse the full richness of angular distributions 
in S — )■ K*{-^ K7r)fi'^fi~ decays. We shall show below that all previously studied observables, 
for example the forward-backward asymmetry A-p-Q, can be easily expressed in terms of our 
new observables. S^""^ and aI""^ are normalized to the CP-averaged dilepton mass distribution 
to reduce both experimental and theoretical uncertainties. Taking the CP average means that 
CP- violating effects in the S'f'^'* are washed out, resulting in a cleaner observable. Taking the 
CP asymmetry, on the other hand, means that any non-standard CP violation can be easily 
identified. 

These CP asymmetries, i.e. Af^\ are expected to be small in the SM, as previously noted 
in Ref. [20]. This is because the only CP- violating phase affecting the decay enters via A„ in 
Eq. (12. ip and is doubly Cabibbo-suppressed. Therefore we are particularly keen to examine 
these asymmetries in the context of CP-violating phases in NP models. 

It should be stressed that out of these 24 observables, two vanish in the SM, namely Sq 
and Aq, which are generated only by scalar operators, and four are related in the limit of 
massless leptons through = 3S2, = —S2 and Af = 3A2, A\ = —A^ (see Sec. 13. 3p . 
Table summarizes the number of independent observables in these limits. 

In addition, even for non-zero lepton mass, only three of the four S^{\ are independent, 
which can be seen as follows. The dilepton mass distribution can be expressed in terms of 
angular coefficients as 

|^ = ^(2/i + /r)-^(2/| + /2). (5.3) 
Therefore, due to the normalization (15. ip . there is the relation 

l{^2S{^SX)-\{:iSf2^S'2) = \. (5.4) 

* Note that our definition of the CP asymmetries differs from Ref. [30] by a factor of |. 
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Consequently, the complete set of 24 independent observables would be given by the twelve 
Af\ eleven S'f^ and the CP -averaged dilepton mass distribution d{V + V)/dq^. However, the 
latter is the only observable for which the normalization of the form factors is relevant, so 
theoretically it is not as clean. 

In our opinion, the quantities s\"''^ and A^f' are the natural starting point for an exper- 
imental analysis. In Ref. [21], a detailed investigation was carried out showing that a full 
angular fit was the preferred way to extract observables. This would involve fitting Eqs. (13. 9p 
and (13. lip to data. From such a fit the and would be found directly, and could be 
combined using Eqs. (15. ip and (15.21) to give the desired quantities. We suggest that a similar 
full angular fit could be carried out for the four- fold spectrum d^iV ± F), so S\"''' and 
would be instantly accessible. Note that, due to Eq. (I3.12p . the CP-averaged decay distri- 
bution (i^(r + f) gives access to S'fl^^-j and Ag^ggg, while the remaining observables can be 
obtained from df'iV — T). 

Alternatively, S\"''' and A^f' can be found by taking asymmetries and/or integrating (i^(r± 
r) over the angles 6i, 6k and 0. Details for the extraction of some of the A^f^ are given in 
Ref. [20], but we stress that all our observables can be determined in a similar manner. To 
illustrate this point, one case not mentioned in Ref. [20] is 5*5, which can be obtained by 
integrating over two angles: 



= — r 



37r/2 



tt/2 



tt/2 



2tt 



3tt/2 



d cos 9 k 



d\r-r) /d(r + r) 



dq"^ d cos 9kc 



dq^ 



(5.5) 



As stated above, we normalize the S^'^^ and a['^^ to the CP-averaged dilepton mass dis- 
tribution in order to reduce the dependence on the form factors. Our approach described in 
Sees. 13.21 and [3.41 makes use of the full form factors for the dominant leading-oder contribution 
and the soft form factors for additional suppressed contributions. Therefore our results are 
largely independent of the relation between the soft form factors and the full form factors. 
However, to further our understanding of these soft form factor relations, we investigate them 
and their dependence in App. [Bl It is found that relations involving ^± are almost inde- 
pendent of g^, whereas those involving have a considerable dependence on due to the 
neglected 1 / mi, terms. Therefore we stress that the transversity amplitudes A^^ of Sec. 13. 2|, 
and all angular observables built from them, should be more or less insensitive to 1/mi, cor- 
rections, i.e. corrections to QCDF, while Aq'^ and all corresponding angular variables will be 
slightly more affected by such corrections. These findings impact on prior work carried out in 
this channel, where the transversity amplitudes were given entirely in terms of the soft form 
factors using QCDF. 



All established observables can be expressed in terms of S^^^ and Af^' . For example, the 
CP asymmetry in the dilepton mass distribution is given by (see Eq. (15. 4p ) 



A 



CP 



d{T-T) /d{T + T) 



dq" 



dq^ 



^(2A? + A5)-1(2A^ + A^). 



(5.6) 



We prefer to define the normalized forward-backward asymmetry as a ratio of CP-averaged 
quantities, to wit 



A 



FB 



d cos Ui 



dq^ d cos 9i 



'^f(r + r) _3 

dq^ 



SI). 



(5.7) 
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The CP average is numerically irrelevant in the SM, but makes the connection to experiment 
more transparent. In addition, this definition is complementary to the forward-backward CP 
asymmetry [79] . 



4 CP 



' dq^ d cos 9i I d(f 



Additional well-established observables are the K* longitudinal and transverse polarization 
fractions F^, Ft, which are usually defined in terms of transversity amplitudes. We prefer to 
directly express them in terms of CP-averaged observables and define 

Fl = -SI Ft = (5.9) 

The well-known relation Ft = 1 — Fl is then a consequence of Eq. (15.41) in the limit of 
vanishing lepton mass. 

In Refs. [IHIET], the transverse asymmetries A^^ were introduced. They can be expressed 
in terms of our observables as 

,(2) _ ^3 



2 



(3)_/ 4 51 + 52 



-2SS(2SS + S, 



Finally, for some observables it is useful to consider their average. We define 

6GeV2 , , / ^6GeV2 



d{T + T) 

aq 

IGeV^ ^ ' / JlGcV^ 



I dq'{lt^ + lt^) I I dq^^-"^^^ (5-11) 



4'"' =/ diil't'-lf] df'^S^. (5.12) 



IGcV 



2 



iGcv2 dq^ 



The reasons for choosing the interval 1 GeV < < 6 GeV are discussed in Sec. |2l 

We proceed in the next section by studying the predictions for s\'^^ and A^f* , keeping in 
mind the sensitivity to hadronic effects. This is carried out first in the SM and later in the 
various NP models described in Sec. HI 



6 Phenomenological Analysis 

We are now in a position to perform a phenomenological analysis of the observables defined 
in Sec. first in the SM, then in a model-independent manner, and finally for specific NP 
scenarios. 
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6.1 Standard Model 



Our predictions for the CP-averaged angular coefficients in the SM are shown in Fig. |2l 
Sf and Si have been omitted since the relations Sf = SS"! and S'f = —S2 (see Sec. 13. 3p are 
fulffiled up to lepton-mass effects, which amount to at most 1%. Si 2 are numerically large 
as expected. 5*4, 6*5, Sq are similar in magnitude, but are particularly interesting as they 
each have a zero in g^. All these predictions are seen to have small uncertainties, as the 
normalization results in a cancellation of hadronic effects. In Tab. |6l we show our predictions 
for the positions of the zeros of 5*4, 5*5 and S'g, denoted by QoiSl""^) from now on. S3 is 
numerically small in the SM since it is approximately proportional to the chirality-flipped 
Wilson coefficient C7, which is suppressed by a factor m^/mf,. 5*7, 5*8 and Sg are small as 
well and have a larger error-band as they arise from the imaginary part of the transversity 
amplitudes. 

The error bands have been obtained by adding various uncertainties in quadrature. We es- 
timate the uncertainty due to the form factors by varying the Borel parameter and continuum 
threshold as discussed in Sec. 12.21 The renormalization-scale uncertainty is found by varying 
fi between 4.0 and 5.6 GeV, where fi is the scale at which the Wilson coefficients, a.^ and 
the MS masses are evaluated. We also include parametric uncertainties which are estimated 
by varying the hadronic parameters as indicated in Tab. [HI the ratio rric/mh between 0.25 
and 0.33, and the CKM angle 7, which is particularly important for the doubly Cabibbo- 
suppressed contribution to the CP asymmetries, between 60° and 80° |^ In addition, we show 
the leading-order prediction as a dashed line. We find that the impact of radiative QCDF 
corrections is moderate for observables like 5*2,3,4,5,6 that, in the SM, are largely independent 
of weak or strong phases, but becomes more prominent for observables built from imaginary 
parts, like 6*7,8,9 and Ai, where the main contribution comes from strong phases induced by 
0{as) corrections in QCDF. 

Some of these Sl"''^ can be directly compared to previous results in the literature. S"! and 
5*2 correspond to the K* longitudinal and transverse polarization fractions and Ft, see 
Eq. (15.91) . and Sq yields the forward-backward asymmetry Ap^, see Eq. (15.71) . In particular, 
QqISq) in Tab. [6] is identical to the zero of the forward-backward asymmetry which has been 
extensively studied in the literature. For completeness, in the last row of Fig. [2] we also show 
the CP averaged dilepton mass distribution d(r + f )/ dq^ and the observables A^^^ and Aip^ 
defined in Ref. [21], see Sec. |5l We find that our results for all these observables compare well 
to those in the literature. However, we note that the peak in the plot of is a manifestation 
of the zero qQ{Si) of 5*4, see Eq. (I5.10p . This division by a near-zero quantity induces a large 
theoretical uncertainty both in the position of the peak and its height. We stress that such 
uncertainties do not arise if the observables S4 and S^ are considered instead of and A^\ 
In fact, as dV /dq^ is a smooth function in the range of q^ considered, none of our observables 
Si and Ai is affected by accidental and delicate cancellations in the denominator. 

As explained in Sec. [5], the CP asymmetries are close to zero in the SM, which is evident 
from Fig. [3l where we show all the A^f* (again except for AJ'*^) and the CP asymmetry in 
the decay distribution, Aq^. As explained above, the shift from LO to NLO is substantial. 
Our results are in good agreement with Ref. [2D], but do not coincide exactly. This can be 
understood by recalling that we use the full LCSR form factors and that our normalization of 

^The discontinuity in some of the error bands just below 6 GeV^ is an unphysical artifact resulting from 
small charm quark masses ~ 1.2 GeV allowed in the estimation of the error. This feature was already observed 
in Ref. [12]. 
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Figure 2: CP-averaged angular coefficients s\"'\ CP-averaged dilepton mass distribution d{T + 

T)/dq^ and transverse asymmetries ^^'^^ in the SM as a function of q'^. The dashed lines are the 
leading-order (LO) contributions, obtained in naive factoriation. The thick solid lines are the full 
next-to-leading order (NLO) predictions from QCD factorization (QCDF), as described in Sec. 12.41 
The blue band defines the total error for the NLO result as described in the text. 
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Figure 3: CP asymmetries A^f^ and Acp in the SM as a function of cp'. The meaning of the curves 
and bands is as in Fig. [2j 
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the soft form factors, especially is different from that used in Ref. [20]. Also our choice of 
quark masses, in particular nic/mi, and m^, as well as the scale /i at which the QCDF hard- 
scattering corrections are evaluated, differs from [20]. We stress that, in view of the smallness 
of the SM values of Ai, these discrepancies become irrelevant once large NP contributions 
start to dominate these observables, as we shall see in the remainder of this section. 

In Tab. [71 we list our predictions for the g^-integrated CP-averaged angular coefficients 
and CP asymmetries as defined in Eqs. (15.111) and fl5.12p . {S2), (5'g) and (Aqp) can be 
directly compared to existing experimental results from BaBar and Belle [6l [91] . 

6.2 Model-independent Considerations 

Before turning to specific NP scenarios, we investigate the model-independent impact of the 
Wilson coefficients on our observables. 

6.2.1 Impact of Wilson Coefficients on Observables 

The impact of NP on the angular observables discussed in our paper is given by the changes 
of the Wilson coefficients of the affected operators. One can group these Wilson coefficients 
into three classes: 

• Dipole coefficients: C7, C7, Cg and Cg. The role of the gluon dipole operators is sub- 
leading in the decay considered. 

• Semileptonic coefficients: Cg, Cg, Cio and C[q. 

• Scalar coefficients: Cs — Cg and Cp — C'p. 

Before entering the discussion of various NP scenarios, it is useful to study the correlation 
between the angular coefficients and the Wilson coefficients. In Tab. [8| we show which ob- 
servables are most affected by a significant change of a given coefficient. In Tab. [9] we show, 
on the other hand, which Wilson coefficients should be altered to produce a large effect in 
specific observables. 
We observe: 

• C7, C7, Cg, Cg, Cio and C[q can induce large effects in many observables, or at least 
in those that do not require the presence of strong phases. To be precise, the Ai are 
mainly induced by imaginary parts of the Wilson coefficients, while the Si are induced 
by their real parts. 

• Only the primed coefficients C7, Cg and C[q can significantly affect the observables 5*3 
and Ag. As can be seen from Eq. f lS.lOp . 6*3 corresponds to the transverse asymmetry 
Aj, and the impact of NP physics contributions to C'-j on this observable has been 
studied for example in Refs. [211 HSl HZ] • 
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Table 6: Predictions for the zero positions q^iSi ) of ^4, ^5 and Sg in the SM. 
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Table 7: Predictions for the integrated CP-averaged angular coefficients {S^""^) (in units of 10^^) and 
the integrated CP asymmetries (^4^^^^) (in units of 10~^) in the SM. Note the different normalization 
of the (j4 ."^ ) with respect to Ref . [20] , see footnote El 



Wilson coefficients 


largest effect in 




cs QC OS QC Q Q Q QS 
'-'l) '-'2' '-'2' "-"S) ^4) *J5) "-"e' 

A, ^, ^, 

BR(S ^ X,7), BR(S ^ 


'-'95 <-^95 '-^10! "-^10 


cs QC OS QC Q Q Q QS 
'-'D '-'d '-'2' '-'2' '-'3) ^4; ^5) ^6' 

Ay, As, Aq, 
BR{B 


Cs — C'g 


CC 


Cp — Cp 


BR(5, ^ fi+fi^) 



Table 8: The Wilson coefficients relevant in i? — t- K*fi^fi and the observables they have the largest 
impact on. 



• The scalar operators affect mainly and the branching ratio for Bg — )■ . This 
impUes interesting correlations between these two observables as discussed in Sec. I6.2.3I 



6.2.2 Model-independent Analysis of ^4, ^5 and 

The zero of the forward-backward asymmetry has been the focus of many experimental and 
theoretical studies (see for example Refs. [TOl [T7]) as it is established as being an observable 
free from hadronic effects and capable of distinguishing between NP scenarios. In Sec. [5] we 
expressed the CP-averaged forward-backward asymmetry in terms of Sq through Eq. (15. 7p . 
so Sq could clearly be studied instead of ApB- In addition, from Fig. El we find there are 
two more observables with such a zero in g^, 6*4 and S5. A study of these three observables 
in a model-independent way could allow us to constrain the NP contributions to the Wilson 
coefficients. 

From Tab. IHwe see that the zero of S4, Qq^SD, is largely sensitive to C7, Cj, Cio and C(q. 
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Table 9: The most interesting anguiar observables in B 
tfiey are most sensitive to. 
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Figure 4: Left: correlation between q^{S/i)-, the position of the zero of ^4, and the NP contribution 
to Cio — C'lQ. Right: correlation between (?o('^6) ^"^^ contribution to Cg. We use the 

branching ratio for B — )• X^^ to constrain the NP contributions to C7 and C^. The green (red) band 
corresponds to a value of BR(i? — )■ ^57) at the upper (lower) end of the experimental 2a range, the 
blue band to SM values for C7, C7. 



This dependence arises only through C7 — C7 and Cio — C[q. We therefore explore how the 
position of the zero in is affected by NP modifications to Cio " CJg and C7. The current 
experimental value of the branching ratio of 5 — provides a constraint on C7 and C7. 
We find a strong dependence of q^i^Si) on Cio ~ C'lQ' its measurement would provide very 
interesting information about these Wilson coefficients. In Fig. |U we show this dependence 
for real values of C7. 

If the NP introduces an imaginary part to C7, the bound from B — )■ Xs'^ is weakened, 
allowing large effects in the zeros. In fact, large values of Im(C7) significantly enhance the 
branching ratio of the decay B — )■ ^^7 and in order to be in agreement with the experimental 
data, large positive contributions to Re(C7) that interfere destructively with C7'^ are required. 
For such values of the Wilson coefficients, the branching ratio of the decay B — )■ Xg^x^ ijT is 
largely enhanced, effectively setting a new upper bound on Re(C7). In the left-hand plot in 
Fig. O we show these combined constraints on C7 in the complex plane. Exactly the large 
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Figure 5: Left: Experimental constraints on the NP contribution to Cy. The blue circles show the 
constraint from the central and ±2(T values of BR(i? — )• Xgj), assuming = 0. The black circle 
corresponds to the 2a bound from BR(i? — t- Xs£^i~), assuming Cfj^^ = 0. The sohd thick and 
the dotted lines have been obtained assuming SM and SMib25% values for Cg, respectively. Right: 
Correlation of the zero in Sq with the NP contribution to Re(C7). The blue, red and green bands 
indicate SM, SM+25% and SM— 25% values for Cg with the associated theoretical uncertainty. The 
vertical dashed lines correspond to the upper and lower bounds on Re(C7) in the absence of an 
imaginary part of Cy. (The corresponding points in the left-hand plot are highlighted by red dots.) 
For an arbitrary imaginary part, the upper bound on Re(C7) is removed, and 0^(3^) can be at or 
below IGeV^. 



positive contributions to KelCj), which are allowed in the presence of phases in Cf^, then 
unambiguously shift the zeros of 5*4, and Sq towards lower values. In the right-hand plot 
in Fig. m we show as an example that the allowed range for q'o(5'g) is greatly enhanced in the 
case of complex Cy. 

This analysis can also be applied to S^, which depends strongly on Cr and Cg. We examine 
the dependence of q'ol'S'e) on NP contributions to Cg and C7. This again is restricted by the 
experimental value of the branching ratio of B ^s7- We find a strong dependence on Cg, 
and for real values of C7 this would be a clean way to determine infomation about a possible 
NP contribution to Cg as seen in Fig HI Again, if NP induces a complex phase of C7, the 
range in QqISq) increases dramatically. 

It is a greater challenge to extract information about the Wilson coefficients from 5*5 due 
to its dependence on C7, Cy, Cg and C[q. However, a measurement of qQ^S^) could provide 
a consistency check with Cio — C[q and Cg determined from 5*4 and 5*6, provided C7, Cy are 
real. In addition, this might allow one to untangle the effects of Cff^ and C'^^ in Fig. HI 



6.2.3 Impact of Scalar Currents 

As mentioned in the introduction, the impact of the scalar and pseudoscalar operators Og^p 
on the angular distribution of 5 — K*{^ K7r)fi'^n~ has been considered before [13], and no 
relevant effects on the observables of interest were found. However, as shown in Sec. 13. 2[ the 
inclusion of lepton-mass effect j^. which were neglected in previous studies, gives rise to an 

^"We stress that we restricted ourselves to muons in our numerical analysis. 
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additional observable in models with scalar currents, which can serve as a precision null-test 
of the SM and, as we will show, in principle allows one to distinguish between different NP 
models. 

To assess the size of the possible effects generated by these operators, we first consider 
the allowed ranges for the Wilson coefficients C^'p. The most stringent constraint on these 
coefficients comes from the measurement of Bg — >■ which is strongly helicity suppressed 

in the SM, with a predicted branching ratio of [5T| [52] 



BR{Bs = (3.37 ± 0.31) x 10' 



(6.1) 



The most recent experimental upper bound still lies, at the 95% confidence level, one order 
of magnitude above the SM [80] : 



BR{Bs fi^fi") < 5.8 X 10" 



(6.2) 



However, in many models, e.g. the MSSM at large tan /3, this branching ratio can be greatly 
enhanced. 

In a generic NP model, the branching ratio is given by 



where 

S -- 



Ami 



m 



Bs 



\S\ 



Ami 



m 



Bs 



\p\ 



m 



Bs 



(Cs-C's), 



P 



m 



Bs 



{Cp-C'p)+m^{C,o-C[o). 



(6.3) 



(6.4) 



2 V - 2 

Considering the experimental bound in Eq. f l6.2p . these formulae imply the approximate 
bounds 



C^l < 0.12 GeV" 



-0.09 GeV"^ <Cp-C'p< 0.15 GeV"\ 



(6.5) 



barring large NP contributions to the Wilson coefficients C^'^ 



10 • 



Now, inspecting the formulae for the angular coefficients, Eqs. fl3.34p - fl3.45p . one can see 
that the only terms in which Cg^ and are not suppressed by the lepton mass enter in the 
angular coefficient Jf . However, due to the small size of the Wilson coefficients themselves, 
see (16.51) . these terms turn out to be numerically irrelevant in general once the bound from 
Bg — 7- fJ^~^fi~ is taken into account. 

Since the pseudoscalar operators do not contribute to any other angular coefficient, this im- 
plies that they are indeed irrelevant in the phenomenological study of -B — >■ K*{^ K7r)fi~^fi~ . 
For the scalar operators, however, the situation is different, because of the new angular co- 
efficient Jg, Eq. f l3.4ip . which is directly proportional to the real part of {Cs — C'g) and thus 
vanishes in the SM. So, although numerically small, this angular coefficient is an appealing 
observable because any measurement of a non-zero value would constitute an unambiguous 
signal of scalar currents at work. 

This is in contrast to the process Bs — )■ where a large enhancement of the branching 

ratio compared to the SM could be caused by both scalar and pseudoscalar currents. In 
addition, the measurement of a non-zero Sq (the CP-averaged counterpart of Jg) would allow 
to determine the sign of Re(C5 — C'^). In fact, by a combined study of Bg — /i and the 
observable Sq, one would be able to constrain the relative sizes of the scalar and pseudoscalar 
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10' (SD 

Figure 6: Correlation between the observable (5g) and the branching ratio of Bg — )■ /U'^/^^- The blue 
band is obtained by assuming NP contributions only to the Wilson coefficient Cs, the black curves 
(where error bars are omitted) by assuming Cp = —Cs- Different values of the phase Arg(C5) are 
indicated. The red and green dots correspond to points in the CMSSM as described in the text. 
The horizontal dashed lines indicate the SM prediction for BR(i?s — t- (j6.ip and the current 

experimental upper bound (|6.2p . 



Wilson coefficients, which can serve to distinguish different models of NP. For example, in 
the MSSM, the ratio of Cs and Cp is 



Cp^ M% 



Cs Ml, 



(6.6) 



to a very good accuracy, a relation which could be tested by a measurement of BR(i?5 — )■ 
/i+/i~) and SI- 

To illustrate this point, we show, in Fig. [6l the correlation between BR(i?<j — > /i^yU~) 
and (5*0 ) (as defined in Eq. (15. lip ). The blue band has been obtained by assuming that NP 
contributions enter only through Cg, i.e. setting Cp/Cs = 0, and varying Cs accordingly; the 
error band takes into account all the sources of error as discussed in Sec. 16.11 

Assuming, in contrast, Cp/Cs = — 1? as would be the case in the MSSM, one obtains 
the black dashed parabola. As an illustration, the predictions for parameter points in the 
constrained MSSM (CMSSM) with large tan/? are indicated as red and green dots. These 
points have been generated by a random scan of the CMSSM parameters in the ranges 

mo < 1 TeV, mi/2 < 1 TeV, (6.7) 

-2mo <Ao< 2mo, 30 <tan/3 < 50, (6.8) 

permitting both signs for the /^-term and discarding points violating existing mass bounds 
or being incompatible with the measurement of BR(i? — ^^7). The green dots correspond 
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to /i > 0, the red ones to < 0. It can be seen that the CMSSM points he on the curve 
corresponding to Cs = —Cp and, in particular for a positive /i parameter, could be clearly dis- 
tinguished from the scenario without pseudoscalar currents, assuming sufficient experimental 
accuracy. 

Since the observable (S'g) probes the real part of {Cs — C^), the correlation gets modified 
if one allows a phase in Cs- More precisely, | (S'g) | gets reduced for a fixed value of BR(_Bs — )■ 
This is illustrated by the black curves corresponding to Cs = —Cp, where both Wilson 
coefficients are now complex, with the respective phase Arg(C5) indicated by the labels on 
the curves. This is precisely what happens in the Flavour Blind MSSM discussed in detail 
in Sec. 16.3.21 In this scenario, the measurement of the correlation between BR(i?s — ?■ 
and (Sq) would thus directly probe the phase of the scalar Wilson coefficient. 

To summarize, while pseudoscalar operators are numerically irrelevant in the decay B — )■ 
K*{—^ K7i)fi'^ fi~ , a study of the angular distribution allows one to probe the scalar sector of a 
theory beyond the SM, in a way that is theoretically clean and complementary to Bg — )■ fJ'~^fi~- 

6.3 Specific New Physics Scenarios 

With the SM predictions for the CP-averaged angular coefficients S^""^ and the CP asymmetries 
A^""^ in hand, we now investigate how these observables change in the NP scenarios discussed 
in Sec. HI 

6.3.1 Minimal Flavour Violation 

In the MFV framework as described in Sec. 14. 3[ no additional CP-violating phases are present 
and NP contributions to the Wilson coefficients of the primed operators can be neglected. This 
implies that all possible effects will arise from real contributions to the Wilson coefficients 
Cy, Cs, Cg and Cio- This implies in turn that the most visible departures from the SM 
predictions will be in the observables 5^ 2, 5*4, S5 and Sq, while the other angular observables 
and in particular the CP asymmetries will essentially be SM likel"1 

Model-independent studies within the MFV framework show that large NP contributions 
to the Wilson coefficients C7, Cg, Cg and Cio are still allowed [81]. In particular, scenarios 
in which the sign of these Wilson coefficients is flipped with respect to the SM cannot yet be 
excluded. 

However, in concrete MFV NP models it is usually difficult to generate large effects in Cg 
and CiQ. For example in the MFV MSSM, NP contributions to Cg and Cio are typically very 
small [76l [82]. Therefore, in this model, the main source of NP effects is Cy whose value can 
be modified substantially by chargino-stop loops. For negligible NP contributions to Cg and 
Cio however, the effects in C7 are strongly constrained by the data on BR(i? — >■ ^^^7) and 
BR{B XsH'^n') and in particular a sign flip in C7 is excluded at the 3cr level [53]. The 
effects in the S\'^^ are then quite limited. In Fig. [3 we show the largest possible effects in S'4, 
5*5 and S^. scenario MFVi (green curves) corresponds to the maximum allowed negative (i.e. 
constructive) NP contribution to C7 (i.e. Cf^) and shifts the zeros of S'4, 5*5 and to larger 
values of q^. Scenario MFVn (red curves), on the other hand, corresponds to the largest 
positive allowed value of C^^ and hence shifts the zeros to smaller values. The separation in 

"'^^ In general the concept of MFV does not exclude effects in the scalar Wilson coefficient Cs which affect 
the observable Sq. However, as shown in Sec. 16.2.31 these effects can be discussed completely independently 
and we will not consider them in this section. 
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Figure 7: The observables S^, S5 and in the SM (blue band) and the MFV MSSM scenarios 
MFVi n described in the text. 
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Figure 8: The correlation between the zeros of S5 and 5*1 and BR(i? — ^ ^s"}) in the MFV 
MSSM. The blue circles correspond to the central SM values, while the green diamonds represent 
scenario MFVj and the red squares scenario MFVn. 



between these two curves corresponds to the range shown in Fig. Hlfor {Cf^f — C'^^) = 
and C^^ = 0, respectively, where the superscript NP denotes the NP contribution to the 
Wilson coefficient. The most relevant input parameters corresponding to the two scenarios 
are collected in Tab. [101 



Scenario 


tan P 


niA 


rrig 






At 




MFVi 


28 


380 


530 


800 


540 


-850 


860 


MFVii 


29 


530 


1000 


880 


660 


880 


750 



Table 10: Most relevant parameters of the two MFV MSSM scenarios discussed in the text, tan (3 
is the ratio of the two Higgs VEVs, uia the mass of the pseudoscalar Higgs, rrig is the gluino mass, 
TTig is a universal soft mass for the left handed squark doublets, rrijj a universal soft mass for the 
right handed up squarks, is the stop trilinear coupling and /i the Higgsino mass parameter. Our 
conventions for the trilinear coupling are such that the left-right mixing entry in the stop mass 
matrix is {m?)LR = —nitiAj: + /i* cot P). All massive parameters are given in GeV. 

It is well known that in the MFV MSSM, the shift in the zero of the forward-backward 
asymmetry in S — )■ Xsfi~^fi~ is highly correlated with a change of the branching ratio of 
B — >■ Xs'~f [551 El]- In Fig. IHlwe show the corresponding correlation between the zeros of S4, 
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Scenario 


tan (3 


rriA 








At 




Arg(/iA-) 


FBMSSMi 


40 


400 


700 


380 


700 


900 


150 


-45° 


FBMSSMn 


40 


400 


700 


380 


700 


900 


150 


50° 


FBMSSMiii 


40 


400 


700 


650 


700 


900 


150 


60° 



Table 11: Most relevant parameters of the three FBMSSM scenarios discussed in the text. Ah 
massive parameters are given in GeV. 

5*5 and Sq and BR(i? — )■ X^^). Any deviation from the lines in the plots would signal the 
presence either of NP contributions to Wilson coefficients other than C-i or of new CP- violating 
phases that lead to complex values of C7. 

6.3.2 Flavour Blind MSSM 

One model with new sources of CP violation is the FBMSSM discussed in Refs. [SHI EOl EB E2] ■ 
This is a MSSM where the CKM matrix is the only source of flavour violation, but additional 
CP-violating, flavour conserving phases are present in the soft sector. Within this framework, 
the majority of non-standard effects arises though complex NP contributions to the Wilson 
coefficient Cj. We discuss two scenarios in which the effects are maximal: scenario FBMSSMi 
is characterized by large negative Im(C7), while scenario FBMSSMn corresponds to a large 
positive Im(C7). The corresponding input parameters are collected in Tab. [HI together with 
those of a third scenario, FBMSSMni, to be considered later. 

Concerning the CP asymmetries, we observe that significant departures from the SM 
predictions can be obtained in AI2, A^, A^, Aj and Ag. The most pronounced effects can 
be seen in Aj and Ag and these are shown in the left and centre plot of Fig. |9l The effects 
here are predominantly due to the large imaginary part of Cj and we note that in this case 
positive values for Aj imply negative ones for A^ and vice versa. This is also displayed in 
the right plot of Fig. |9l where we show the almost perfect correlation between the integrated 
asymmetries {Aj) and {As). Any deviation from the line shown in this plot would signal the 
presence of additional imaginary parts in either Cj or Cg^ and Cj'J. 

In the CP-averaged angular coefficients we find significant departures from the SM in 812, 
5*4, 5*5, Sq and also in Sq, while effects in 5*3, 5*7, 5*8 and Sg can hardly be distinguished from the 
SM. Although in the FBMSSM the BR(i?<. — )■ fi~^fj^~) can be close to its experimental upper 
bound, the effects in are smaller than the maximal effects found in the model-independent 
discussion of Sec. I6.2.3[ because the large imaginary part in C7 implies a large phase for the 
relevant Wilson coefficient Cs- Concerning Sl'2, we find that while jS"! 2I is enhanced, jS'^jl is 
suppressed with respect to the SM results. For S^, and the forward-backward asymmetry 
Sg we find significant shifts in their zero towards values of lower than the SM prediction 
or we even find no zero at all. These effects are shown in Fig. [TU] and are much larger than 
those possible in the MFV MSSM (see Fig. [7]). The reason for these large shifts are the large 
values of Im(C7) in the scenarios considered, as discussed in Sec. 16.2.21 

One finds that the strict correlation between the zeros and BR(i? — )■ Xg'-f) is lost in 
the FBMSSM. This is shown in the upper plots of Fig. [TT] However, as the additional 
contributions to 6 — )■ S7 from the imaginary part of C7 can only enhance the branching ratio, 
one still finds an upper bound on the zeros for a given value of BR(i? — >■ Xg'j). In addition, in 
the lower plots of Fig. [11] we show the zeros gQ(S'4), qoiS^) and Q'ol'S'e) against the integrated 
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Figure 9: Left and centre plot: CP asymmetries Aj and in the SM (blue band) and three 
FBMSSM scenarios as described in the text. Right plot: correlation between the integrated asym- 
metries (^7) and (^g) in the FBMSSM. Blue circle: SM, green diamond: FBMSSMi, red square: 
FBMSSMii , orange triangle: FBMSSMm. 
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Figure 10: The observables 5*4, ^5 and in the SM (blue band) and the three FBMSSM scenarios 
FBMSSMi h m. 



asymmetry (Ay). One observes that large effects in (Ay) are correlated with large shifts in 
the zeros towards lower values. 

In order to identify signs in the CP asymmetries which are favoured in this model one 
must include additional observables in the analysis. To this end we also investigate the direct 
CP asymmetry in the 6 — >■ S7 decay Acp(fe — 57), the electric dipole moments of the electron 
and the neutron df. and (i„ and the mixing induced CP asymmetry S^Ks- We recall that 
in [62] striking correlations between these observables have been found. In particular, the 
desire to explain the anomaly observed in S^^Ks through the presence of flavour conserving 
but CP- violating phases implied a positive Acp(& — 57), by an order of magnitude larger 
than its SM tiny value and dg, at least as large as 10"^^ ecm. 

The left plot of Fig. [T^ shows the correlation between (Ay) and S^Ks- We find that a value 
of S(i,Ks — 0.44, as indicated by the present data [85], implies a negative value for {Aj) in the 
range [—0.2, —0.05] and then also a positive value for (As) in the range [0.03, 0.11]. In addition 
to the two scenarios discussed above, we have chosen also a third scenario, FEMSSMm, 
indicated as orange triangle in the plots of Figs. [9l [11] and [121 that gives S^Ks close to the 
experimental value. This scenario is shown in Figs. M and [TD] as the orange bands and we find 
that while one still can get almost maximal effects in (Ay) and (Ag) the effects in 6*4, 5*5 and 
5*5 are much less pronounced. 



38 



2.1 r 

2 2.0 

X 

^ 1.9 



1.; 



1.0 



0.3 

0.2 — 
0.1 — 
(A^) 0.0- 
-0.1 — 
-0.2—; 



ql{S^) (GeV2) 



qliS,) (GeV2) 



2.5 



o 


2.0 


X 






1.9 


T 


1.8 










m 


1.7 



'1o 



0.3 

0.2 — 
0.1 — 

{A^) 0.0- 
-0.1 — 
-0.2 — 



w 



1.5 2.0 2. 

gg(55) (GeV^) 



1.5 2.0 2.5 
qi(S,) (GeV2) 



2.1 



^ 2.0 
X 



1.9 



1.8 



1.7 



1.5 



1 

0.3 
0.2 
0.1 

(Ay) 0.0 
-0.1 
-0.2 
-0.3 



1 



3 4 
2s 



qiiSl) (GeV^) 



9^(5^) (GeV2) 



Figure 11: Correlation between the zeros of ^4, ^5 and 5g witli tiie 6 — 57 branching ratio (upper 
plots) and with the integrated asymmetry (^7) (lower plots) in the FBMSSM. The blue circles 
correspond to the SM predictions. The orange triangles correspond to a FBMSSM scenario that 
gives S^Ks close to the central experimental value ~ 0.44. 



In the centre plot of Fig. [T2] we report the correlation between {Aj) and Acp(fe — )■ S7). 
One observes that negative values for (^17) imply positive values for v4cp(& — ?■ 57) that can 
reach values up to (5 — 6)%. 

Finally, the right plot of Fig. [12] shows the correlation between {Aj) and the EDM of 
the electron, de in the FBMSSM. We find that large values for (^17) necessarily require large 
values for the electron EDM close to the current upper bound of 1.6 x 10~^^ ecm |86j . 

6.3.3 LHT 

We analyse the angular observables within the LHT by means of a global parameter scan 
taking into account all relevant constraints from other flavour observables. As already antic- 
ipated in Ref . |68] , most NP effects in the observables considered here are found to be small. 



Scenario 
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XL 




m]j 




ad 


ad 


od 




^13 


^12 


LHTi 


1000 


0.5 


565 


1000 


770 


1.60 


2.50 


1.35 


5.70 


4.20 


5.80 


LHTn 


1000 


0.5 


1000 


375 


425 


1.50 


1.00 


4.75 


4.25 


0.60 


2.85 



Table 12: Parameters of the LHT scenarios LHTt tt: Of- and Sf^ are the parameters of the CKM-like 
unitary mixing matrix for the mirror d quarks, rn^jj are the masses of the mirror quarks, / is the 
high energy scale and xl the mixing parameter of the SM top and the T-even top partner. 
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Figure 12: (Aj) vs. S^^Ks (left plot), (Aj) vs. A^^^p (centre plot) and {Aj) vs. de (right plot) in 
the FBMSSM. The blue circles indicate the SM values, while the green diamonds, red squares and 
orange triangles correspond to the scenarios FBMSSMi, FBMSSMn and FBMSSMni, respectively. 



Ay 




(Ay) 



0.010 



0.000 




(GeV^) 



(GeV^) 



-0.006-0.004-0.002 0.000 0.002 
(As) 



Figure 13: Left and centre plot: CP asymmetries Aj and A^ in the SM (blue band) and the LHT 
scenarios LHTiji. Right plot: Correlation between the integrated asymmetries (^47) and {A^) in 
the LHT. The blue circle represents the SM, the green diamond scenario LHTj and the red square 
scenario LHTn. 



In particular Sq, the forward-backward asymmetry, turns out to be very close to the SM. The 
same applies to all other CP-averaged angular coefficients and most CP asymmetries. The 
largest effects relative to the SM are found in Af and as in the SM their absolute values 
are at most 6 x 10"'^ and 5 x 10~^, respectively. We consider two scenarios, LHTi and LHTn, 
with input parameters as given in Tab. [121 In the left and centre plot of Fig. [13] we show 
the corresponding asymmetries A7 and Ag as functions of q^. The blue curves represent the 
SM. The green curves labelled LHTi correspond to a LHT parameter point that gives the 
largest negative NP contribution to Im(C9) and Im(Cio), while the LHTn curves (red) give 
the largest positive contribution. Enhancement of both asymmetries by a factor of three is 
possible for low values of with visible but smaller effects for larger values of g^. 

Still, these significant enhancements are one order of magnitude smaller than those found 
in the FBMSSM. The reason why much larger effects in Aj and A^ are possible in the latter 
model is that large NP contributions to the imaginary part of C7 are allowed, comparable in 
magnitude to the SM contribution. In the LHT model NP contributions to C7 are found to 
be very small [68]. As the effects in Ai and As are therefore dominantly created by Im(C9) 
and Im(Cio), the correlation between the integrated asymmetries (^47) and (As) is completely 
different than that found in the FBMSSM (see the right-hand side plots in Figs. 19) and [T3l) . 

As a side comment, in our numerical analysis we have used the formulae of Ref. [M] 
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Arg(<5rf)|'2'^ 


GMSSMi 
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520 


500 


400 
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380 


800 


750 


470 


0.01 


-135° 


GMSSMn 
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740 


1000 


460 


1000 


390 


1500 


440 


200 


0.03 


60° 



Table 13: Most relevant parameters of the two general MSSM scenarios with large Cj as discussed 
in the text, a universal soft mass for the right handed down squarks, ^^(j) universal trilinear 

couplings for the up (down) squarks and {Sd)32 the left-right mass insertion that generates large 
effects in Cy. Our conventions for the trilinear coupling are such that the left-right mixing entry 
in the sbottom mass matrix is {m?')LR = — + /U*tan/3). All massive parameters are given in 
GeV. 

modified by the additional term found in Ref. [70] which remove the UV cutoff dependent 
terms in Cg and Ciq. This modification decreases the two asymmetries by roughly a factor of 
2 to 3. Whether this is the final result for the LHT model remains to be seen as the structure 
of the full heavy-fermion sector in the LHT model is rather involved and a complete analysis 
is still lacking. 

6.3.4 General MSSM 

Due to the huge number of free parameters in the general MSSM, a comprehensive analysis 
of this general framework is challenging. As a first step we therefore restrict ourselves to 
a framework in which NP effects are created dominantly by complex contributions to the 
Wilson coefficient Cj. Such a situation can easily be achieved in the general MSSM if one 
introduces flavour violating terms only in the left-right sector of the down squark mass. In 
particular, a (Sd)^^ mass insertion will mostly create contributions to Cj by means of down 
squark - gluino loops, while at the same time leaving the other relevant Wilson coefficients 
SM like. 

Fig. [H] shows possible effects in S4, 6*5 and Sq that arise in this framework due to the real 
part of C'j, for two example scenarios, GMSSMi and GMSSMn. In Tab. [T31 we collect the 
corresponding input parameters. Compared to the framework of the FBMSSM (see Fig. [TU]) . 
the shift in the zeros of these observables show a completely different pattern. While the zero 
of Sq remains SM like, a positive shift in qQiS^) implies a negative shift in q^iS^) and vice 
versa. 

Large imaginary parts of C7 lead to sizeable effects in the asymmetries A7 and As, but 
again the pattern of these effects is different to that in the FBMSSM seen in Fig. |9l As 
shown in Fig. [151 ^ positive (negative) Ay implies also a positive (negative) A^. In particular 
the correlation plot in the right panel of Fig. [15] is completely orthogonal to the one in the 
FBMSSM (see Fig. [9]) and thus a clear distinction between these two frameworks is possible. 

In addition a large complex C7 also leads to large non-standard effects in the observables 
5*3 and Ag as shown in Fig. [TBI In fact, as already mentioned in Sec. 16.2. 1^ effects in and Ag 
are characteristic for scenarios with large NP contributions to the primed Wilson coefficients. 
The large effects in S3 are driven by the real part of Cj and directly correspond to the large 
effects in the transverse asymmetry Arp that have been analysed in |l5l [13 [21]. Having 
analysed possible effects in a particular non-minimal flavour violating MSSM framework we 
flnally mention also the case of the general MSSM with generic flavour violating soft terms. 
Instead of presenting an exhaustive discussion of this framework, we concentrate on two 
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Figure 14: The observables 5*4, S5 and 5| in the SM (blue band) and two GMSSM scenarios with 
large complex contributions to C7 as described in the text. 
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Figure 15: Left and centre plot: CP asymmetries A-j and As in the SM (blue band) and two 
general MSSM scenarios with large complex contributions to C7. Right plot: Correlation between 
the integrated asymmetries (^7) and {Aq) in the framework of a general MSSM with large complex 
Cy. The blue circle corresponds to the central SM value, while the green diamond represents scenario 
GMSSMi and the red square scenario GMSSMn. 
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Figure 16: The observables ^3 and in the SM (blue band) and the two GMSSM scenarios 
GMSSMiji with large complex contributions to C7 as described in the text. 
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Figure 17: Several observables in the SM (blue band) and two selected GMSSM scenarios that show 
large non-standard behaviour. See text for details. 
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specific scenarios that show effects that go beyond those discussed in the above. 

Among the Wilson coefficients that are relevant in the decay B — )■ K*fi~^fi~ the ones that 
are most sensitive to NP effects arising from ffavour violating down squark masses are C7 
and Cj. In the plots of Fig. [T7]we show a scenario GMSSMiv that corresponds to large NP 
contributions to both C7 and Cj. In contrast to the scenario with NP effects dominantly in 
Cj discussed above, one observes e.g. sizeable effects in the zeros of and while the zero 
in 5*4 is much less affected. 

One possibility to generate large effects in the Wilson coefficient Ciq in a supersymmetric 
framework is through fiavour violating entries in the left-right part of the up squark mass 
[761 1791 182] . Scenario GMSSMm in Fig. [T71 corresponds exactly to such a scenario where in 
addition to large complex NP contributions to C7 and C'^, Ciq also receives sizeable complex 
corrections through a (5^)32^ mass insertion. These curves show again a qualitatively different 
behaviour in various observables. For example large effects in and Ag can be observed, 
that however do not show a zero in contrast to the red curves discussed above. 

7 Summary and Conclusions 

In this paper we have analysed all angular observables in the rare decay B — )■ K*{^ 
K7r)^~^fi~. They can be measured at the LHC and later at an upgraded Belle and a Super-B 
facility. These angular observables can be expressed in terms of CP-conserving and CP- 
violating quantities and offer new important tests of the SM and its extentions. To this end 
we have improved on previous studies in a number of ways that have been listed in Sec. [1] 

Having identified angular observables with small to moderate dependence on hadronic 
quantities and large impact of NP we have analysed these observables first within the SM 
and subsequently within a number of its extentions like models with MFV, the flavour-blind 
MSSM with new flavour conserving, but CP-violating phases, the LHT model and also within 
a general MSSM with generic flavour violating soft terms. 

The main messages from this study are as follows: 

• The most promising and complete set of observables in this channel are our Sl"^ and 
Af"^ deflned in Sec. |5l 

• Our predictions for the CP-averaged angular coefficients S^""^ in the SM are shown in 
Fig. 121 Some of these S^""^ are found to be large. 

• On the other hand, as evident from Fig. |3l the CP asymmetries Af^^ are close to zero 
in the SM. 

• Our model independent study shows that pseudoscalar operators are numerically irrele- 
vant in the decay B — )■ K*fi~^fi~. On the other hand a study of the angular distributions 
allows, in a way which is theoretically clean and complementary to Bg — >■ /^"""/i", to probe 
the scalar sector of a theory beyond the SM. 

• As one expects, the A^""^ are SM like in MFV models. On the other hand, some of the 
Sl^°'\ in particular 6*4, S5 and can show deviations from the SM as seen in Figs. [7] 
and [HI where results in the MFV MSSM are shown. 
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• Probably the most interesting results are found in the FBMSSM, in which several 
and A^""^ differ significantly, even by orders of magnitude from the SM results, and there 
exists a number of striking correlations among the observables discussed here and also 
correlations between (and Aq) and Acp{b — 57) and S^Ks- All these results are 
shown in Figs. |9] to [T2j 

• In the LHT model only the CP asymmetries Aj and As differ significantly from the SM 
predictions, but these enhancements are smaller than found in the FBMSSM model. 
This different pattern of effects could easily distinguish these two models. 

• As expected, in a general MSSM, the very large space of parameters does not allow 
for clear-cut conclusions. Almost all observables considered in the present paper can 
significantly differ from the SM results and the pattern of deviations can differ from 
those found in the FBMSSM and LHT models. This is illustrated in Figs. [H] to [T71 
This should allow these three models to be distinguished from each other. 

Clearly, it will be very exciting to monitor the upcoming LHC, Belle upgrade and even- 
tually Super-B factory in this and in the next decade to see whether the angular observables 
discussed in our paper will give a hint for any of the extensions of the SM. 
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Appendices 

A Kinematics of Four-Body Decays 

In this appendix we collect some relevant results concerning the kinematics of the four-body 
decay X — >■ y(— )■ a6)Z(— )■ cd). For an excellent discussion see Ref. [87]. The translation of 
the general results to the notation appropriate for the decay B — )■ >■ K7r)fi~^fi~ is given 
in Tab. \M 

The four-body phase space is 4 x 4-dimensional. The on-shell conditions of the final states 
reduce this number to 4 x 3. Moreover, four-momentum conservation eliminates further 4 
degrees of freedom. Eventually, exploiting isotropic symmetry, one can fix the three Euler 
angles and ends up with 5 physical degrees of freedom. It is customary and convenient to 
express them by the following set of variables, introduced first in Ref. [SB] for the decay 
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m^^, the effective mass squared of tlie ab system, rria + rrib < rriah < rnx — rric 



• m^^, the effective mass squared of the cd system, nic + md < rricd < rnx — ^a —^b', note 
that ruab + rricd < rrix] 

• 6y, the angle of the particle a in the cm. system of the particles a and b with respect 
to the direction of flight of (a, b) in the X rest system, < 6'y < vr; 

• 6z, the angle of the particle c in the cm. system of the particles c and d with respect 
to the direction of flight of (c, d) in the X rest system, < 6*^ < vr; 

• 9x, the angle between the plane formed by the decay products (a, 6) in the X rest 
system and the corresponding plane of (c, d), —n < 6x < Ti- 
lt is convenient to combine the four-momenta Pi, i = a,b, c, d, of the final-state particles into 
the following symmetric and antisymmetric momenta: 



Pab =Pa+Pb, 
Pcd = Pc+Pd: 

The diparticle masses are then given by 



Qab =Pa-Pb, 
Qcd = Pc-Pd- 



p2 
^ab 



m 



ab ' 



p2 

^cd 



m 



cd ' 



Now we define the angles 9x, Oy and 9z hy 

Qab ' Pcd 



COS 6 



Y 



ab 



cd 



COS 9; 



Q 



cd ' J^ab 



Q 



cd 



ab 



(A.l) 



sin 9 



Pab X Qafe ) X Pcd X Qcd 



X 



Pab X Q 



ab 



Pcd X Q 



cd 



(A.2) 



It is important to note that the three-vectors in the above definition must be evaluated in the 
respective rest frames of the particles X, Y and Z. 

With the above definitions, it is straightforward to express all remaining invariant products 
of the four-vectors Pab, Pcd, Qab and Qcd in terms of the five variables 9x, 9y, 9z, rrilf^ and 



m 



2 . 

cd' 



PabQab 
PabPcd = P 



m„ — m 



b 5 



PabQcd 

QabQab 
QabQcd 



1 1 

Trie ~ ""^^ 



m 



■p 



cd 



2{ml + ml) 



m 



aocd cos 9z 



cd 



[{ml - ml){ml - ml)p 
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general 


B — i- Knu^fi 


[ao) 


[71 J\) 


(rrl\ 

[ca) 


( n + n-\ 




rriK* 








9 


UY 


Ok* 






"ah 












p 


{k-q) 



Table A: Translation table between the variables of the general four-body decay X — t- Y{-^ ab)Z{-^ 
cd) discussed in App. [Aland the decay B — )• )• KTT)fi~^ . Note that ruab is fixed to uik*, which 
leads to an allowed range of q'^ of (2m^)^ < < {niB — rriK*)^ ■ (3, (3^ and A are defined in Sec. O 



+2aaab{ml — m^) cosOy 
+2aacd{ml - ml) cos 9z 
+4:aab<JcdP cos By cos dz 
+4:aab<Jcdmabmcd sin 6y sin 6z cos 6x] , 

^o^MsPabQibPldQid = -1^^^^ sin Oy sin Oz sin Ox • 

fncd'Triab 

The remaining invariants can be obtained from the above by [ah] o [cd) and 9z -^^ 9y. 
Here we use 

-'- / 2 2 2 \ 

P = 1^{t^X -^ab-^cd). 



Pab = \{rrili, - ml - ml) , 



^ab = \/plb-mlml 



and corresponding expressions with (ab) — t- (cd). Tab. |A] provides the translation of the 
general notations introduced above to the special case B — K*{^ K7i)fj,~^ fj,~ . 
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B Form Factors and QCD Factorization 



As we have seen in Sec. 12.31 the QCDF formulas for B — )■ K*fi^fi~ are expressed in terms of 
two soft form factors C,j_,\\, defined in Eqs. f l2.25p - fl2.26p . rather than the seven form factors of 
Sec. 12. 2[ which imphes certain relations between them, valid in the large energy limit. These 
relations are the topic of this appendix and we shall demonstrate that they are indeed fulfilled 
by light-cone sum rules (LCSRs) - both analytically and numerically. 

Exploiting Eq. f l2.28p for all 7 form factors, one can establish a number of relations which 
are expected to be valid for K* energies E ~ 0{mb) [TTj : 



2EmB 



TUB 



V{q^) 2E V{q^) niB + tuk- 



{tub - 
Air 



rriK* 



In 



nir 



-L 



+ 



agCp rriB AE± 



An AE Cx(g^ 



(B.l) 
(B.2) 



{ruB + mK-)/{2E) Ai{q^) - (m^ - mK*)/mB ^2(9^ 



An 



[-2 + 2L] 



agCp rnB{mB — 2E) 



AFii 



{mB/2E)T2{q^)-T^{q^ 
{mK*/E)A,{q^) 



1 + 



An 



In 



nit 



An 



2 + AL 



{2Ey 



fmBV' 



AEi 



(B.3) 



/i" J An \2EJ {E/mK*)i\\{q^y 

(B.4) 

respectively. The abbreviations on the right hand side are defined as 
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^2 



2E , 2E 
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2E TUB 



AFi 



B \b J 1 



3m 



u 



AFii 



3m 



B 



X 



B 



du 



u 



(B.5) 



in terms of the twist-2 K* DAs and the first inverse moment of the B meson DA, 1/\b- 
E, the K*^s energy, is related to q"^ by 2m^£' = m^ — (terms in m|-, are neglected). 
Numerical values of these parameters are given in Tab. [31 

One requirement for the LCSRs is obviously that they fulfill the above relations, to the 
required accuracy, i.e. in the large energy limit, and that the leading contributions to ^^di) 
come from transversal (longitudinal) DAs. The above relations imply in particular that 
both Ai and A2 must to leading order be given in terms of transversal DAs, while these 
contributions must cancel in the combination (m^ + mK*)/{2E) Ai — {tub — mK*)/fnB A2. 
As this is an important point, we write down the tree-level contributions to these sum rules, 
and that for V, explicitly {u = 1 — u, uq = (m^ — g^)/(so — q"^))'- 



Aoiq') 



^bJb 



du e 



Uo 



2 - 2 1 
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9v{u) 
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1 + u — 

du 



<^{u) 
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AS 



+ •••], 



rrti—U' 



du e "Ji^' 
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du e 



2 - 2 1 
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uA/2 — 
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(B.7) 



— uq 
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V{q' 



TUB + rriK* rrifje'^B/'^' 



mlfB 



2 - 2 1 



' due — r/jt.</)_L(M) 



2(m^ — 



du 



(B.9) 



The dots stand for terms with higher powers of rriK* ■ All DAs multiplying f^^} are longitudinal 
(transversal). Note that, as mentioned in Sec. \2.2\ the suppression factor of higher-twist terms 
is given by rnhniK* / {fril — g^), which indicates the break-down of the expansion for — )■ m^. 
Comparing the above formulas with (IB.3P one finds that for g^ = (i? = 7725/2) the terms in 
(j)± cancel exactly in the combination of Ai and A2 and (1B.3P is reproduced up to terms in 
higher powers of mx* (which are neglected in the large energy limit). For g^ 7^ 0, however, 
there appears to be a non- vanishing term in /^,: 



{niB + mK*)/{2E) Ai - {tub -mK*)/mB A^ 



rrih e 



m%fB 



du e 



2-2 



4>±{u) q^{ml - g^ - u{m% - q^)) 



K* 



u 



2mB{rn?^ — q^] 



u 



. (B.IO) 



At this point it is useful to recall that, ideally, if the continuum model was perfect and all 
terms in the twist expansion were known, the LCSR would be independent of sq and M^. In 
reality, LCSRs are non-trivial functions of both parameters and the form factors are extracted 
within an interval of where the dependence on that parameter is small, i.e. the sum rules 
are evaluated near an extremum in M^, if such an extremum exists. In Ref. [89] we have 
argued that the central value of form factors should be evaluated exactly at the extremum. 
This implies that for both an ideal LCSR with completely known hadronic spectral density 
p{u) and a realistic one with p{u) known to a certain accuracy in the light-cone expansion 
one requires 



d 



dW^ 



2-2 



du p{u) = 



1 2 2-2 9 

i )• / dUe'M^ "Af^ 



"0 



B 



mf 



uq^ 



u 



p(m) = o. 



(B.ll) 



As the integral must vanish, one finds that p{u) effectively equals (m^ — q^)/ ~ q^) p(^) / ^ 
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Figure A: Form factor ratios (jB.ip to ()B.4p . l.h.s. divided by r.h.s., calculated from LCSRs with 
central values of input parameters. In the perfect heavy quark limit, all ratios equal 1. Left: form 
factor ratios based on ^_|_. Solid curve: (jB.ip . long dashes: ()B.2p with TijV ^ short dashes: ()B.2p 
with TijV . Right: form factor ratios based on ^y. Solid curve: dES]), long dashes: (lR4l) . as IS 
evaluated at the scale /i^ = — = 0{mh). 



at and near the extremum. Hence the r.h.s. of ( IB.lOp vanishes at the minimum in if 
the Lh.s. is treated as LCSR in its own right and not as sum of two different sum rules 
(with possibly different extrema in and different optimal values in sq). This is also the 
reason why, in order to obtain values for Ao{0) and ^\\, it is not useful to add the numbers 
for the individual form factors Ai and A2 obtained from LCSRs in Ref. [S^- Instead, the 
above combination of Ai and A2 has to be evaluated anew. We would like to add that the 
cancellation of the contribution of (j)± to ^|[ can also be made manifest by including the factor 
1/E into the dispersion relation, using the above result that factors wm^ — ml + uq'^ return 
zero under the integral if it is evaluated at the minimum in M^. One then can make the 
replacement 

1 1 ml — uq^ 

2E ms ml — q^ 

upon which the contribution in 0^ cancels explicitly so that 

in the LCSR method. Again, the dots denote terms in higher powers of mx*- Similar analyses 
can be done for the other form factors ratios, for instance Ai/V, Eq. ( IB.ip . Making again 



use of the fact that an additional factor u under the integral equals {ml — q^) / {m\ — g^), one 
finds from f lB.7p and (]B.9P that f IB.ip is fulfilled to twist-2 accuracy, but also that the ratio 



deviates from the r.h.s. of f IB.ip at twist-3 level, i.e. at 0{mx* /mh). 

Turning to the 0{as) corrections in (IB.ip to ( ]B.4p . we can confirm that to twist-2 accuracy 
(IB.ip does not receive any such corrections: they cancel exactly between numerator and 



denominator. Reproducing the 0{as) corrections in flB.2p to f lB.4p is less trivial and requires 



to explicitly perform the limit mt — )■ 00 of the integral over u. We refrain from doing this 
analysis explicitly, but refer to Ref. [90] where it was shown, for i? — t- vr transitions, that 
LCSRs fulfill the SCET relations also at 0{as)- 

In Fig. \M we plot the form factor ratios as functions of for central values of the input 
parameters, separately for ratios based on ^_|_ and those based on .^y. All ratios include 
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the 0{as) corrections calculated in QCDF. We find that, overall, the QCDF predictions are 
fulfilled by the full QCDF form factors from LCSRs at the level of 10% or better. Nonetheless 
there is a considerable dependence of the ratios on g^. This is due to l/rrih corrections which 
are neglected in QCDF, but turn out to be, numerically, larger than the factorizable 0{as) 
corrections calculable in QCDF. As a consequence, we expect that the transversity amplitudes 
^±{\\)L,R of Sec. 13. 2[ and all angular observables built from them, should be rather insensitive 
to l/rrih corrections, i.e. corrections to QCDF, while Aol,_r, A, and all corresponding 
angular variables will be affected by such corrections to a larger degree. 
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